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Triaxial  erosion  tests  were  used  to  study  how  the  erosion  rate  in  typical 
embankment  dam  core  materials  is  influenced  by  soil  compaction  and  by  the 
chemistry  of  the  eroding  fluid.  The  tests  were  performed  in  a  specially  modi¬ 
fied  triaxial  cell  which  was  able  to  apply  confining  pressure  to  compacted 
embankment  dam  materials  and  to  direct  the  flow  of  the  eroding  fluid  through  a 
preformed  erosion  channel  to  model  the  conditions  that  exist  in  the  core  of  an 
embankment  dam.  In  addition,  tests  were  performed  to  develop  a  method  for 
testing  the  effectiveness  ol  filters  in  preventing  erosion  from  (Continued) 
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20.  ABSTRACT  (Continued) . 

.cracked  embankment  dam  core  materials.— 

Five  typical  embankment  dam  core  materials  ranging  from  silty  to  clayey 
soils  were  studied.  The  results  of  the  investigation  may  be  summarized  as 
follows: 

a.  Varying  the  compacted  density  had  little  effect  on  erosion  for  silt 
materials  compacted  at  the  optimum  molding  water  content.  For  clay 
materials,  the  erosion  rate  increased  as  the  density  decreased  from 
95  to  90  percent  of  standard  Proctor  maximum  density. 

b.  Varying  the  molding  water  content  had  a  significant  effect  on  the 
erodibility  of  silt  core  material,  but  this  effect  was  less  important 
for  clay  core  material. 

£.  The  minimum  erosion  occurred  when  soil  was  compacted  at  or  slightly 
above  the  optimum  compaction  water  content.  Significantly  higher 
erosion  occurred  in  specimens  compacted  dry  of  optimum. 

d.  Erosion  rates  increased  slightly  with  a  decrease  in  eroding  fluid 
ionic  concentration.  This  increase  was  more  significant  for  silt 
core  materials  than  for  clay  core  materials. 

e.  Tests  performed  to  study  the  effectiveness  of  using  standard  con¬ 
crete  sand  (ASTM  C  33)  as  a  filter  material  for  the  downstream 
protection  of  cracked  core  materials  showed  that  the  filter  sand 
was  able  to  stop  the  erosion  of  even  a  highly  erosive  silt  material 
under  the  action  of  hydraulic  heads  of  up  to  40  m,  which  was  equiva¬ 
lent  to  a  hydraulic  gradient  of  350. 

The  soil  erosion  rate  asymptotically  approached  a  material  and  erod¬ 
ing  fluid  dependent  constant  defined  as  the  Maximum  Erosion  Rate. 

This  constant  facilitates  the  incorporation  of  laboratory  erosion 
test  data  into  embankment  dam  design  criteria. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

By 

To  Obtain 

cubic  feet 

0.02831685 

cubic  metres 

cubic  inches 

16.38706 

cubic  centimetres 

feet 

0.3048 

metres 

gallons  (U.  S.  liquid) 

3.785412 

cubic  decimetres 

inches 

2.54 

centimetres 

pounds  (force) 

4.448222 

newtons 

pounds  (force)  per  square 

foot 

47.88026 

pascals 

pounds  (force)  per  square 

inch 

6894.757 

pascals 

pounds  (mass) 

0.4535924 

kilograms 

square  feet 

0.09290304 

square  metres 

square  inches 

6.4516 

square  centimetrei 

tons  (2000  pounds,  mass) 

907.1847 

kilograms 
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EVALUATION  OF  THE  EROSION  POTENTIAL  OF  EMBANKMENT 


CORE  MATERIALS  USING  THE  LABORATORY  TRIAXIAL 


EROSION  TEST  PROCEDURE 


PART  I:  INTRODUCTION 


Background 

1.  The  erosion  of  impervious  core  materials  is  an  important  problem  in 
Geotechnical  Engineering.  One  early  description  of  a  sudden  and  catastrophic 
failure  due  to  cracking  and  erosion  occurred  in  1864,  when  nearly  a  quarter 
of  the  Dale  Dyke  Dam  in  England  was  breached.  The  waters  flooded  parts  of 
the  city  of  Sheffield  to  a  depth  of  three  meters  and  killed  250  people  in 
the  most  devastating  dam  failure  ever  to  have  occurred  in  Great  Britain 
(Smith,  1971).  A  recent  reminder  of  this  problem  was  the  failure  of  the 
Teton  Dam  in  Idaho  which  caused  the  loss  of  11  lives,  the  loss  of  an  85 
million  dollar  project  and  the  lodging  of  over  400  million  dollara  in  claims 
for  loss  of  private  property  (U.  S.  Bureau  of  Interior,  1977). 


2.  Hsu  (1981)  in  Table  1  provides  a  list  of  49  more  dams  which  have 
experienced  erosion  problems  in  the  last  100  years;  in  33  of  those  dams, 
erosion  occured  through  the  embankment  while  in  the  remaining  dams  ero¬ 
sion  occured  through  the  foundation.  Fortunately,  for  most  of  these 
dams,  there  was  adequate  warning  and  catastrophic  damage  was  avoided  by 
lowering  the  reservoir  water  level  and  by  performing  remedial  repairs. 

Current  State  of  Knowledge 

3.  Available  field  data:  In  the  last  decade,  an  increasing  number 
of  case  histories  have  been  presented  in  the  literature  describing  em¬ 
bankment  dam  seepage  and  erosion  (Bertram,  1967;  Sherard,  1973;  USCOLD, 
1975).  Table  2  summarizes  these  and  other  field  observations  of  dams 
which  have  experienced  erosion  problems  and  provides  useful  quantitative 
design  values. 

4.  As  may  be  seen  in  Table  2,  measurements  of  the  volume  of  leakage 
flow  as  a  function  of  time  are  usually  available.  In  many  cases,  erosion 
problems  are  first  detected  by  noticing  rapid  increases  in  leakage  flows 
and/or  the  appearance  of  longitudinal  cracks  or  crater  like  depressions 
at  the  crest  of  a  dam.  In  most  cases,  these  observations  provided  ade¬ 
quate  warning  time  to  perform  remedial  work  to  prevent  breaching  of  the 
dam. 

5.  However,  there  is  almost  a  complete  lack  of  data  describing  the 
rate  at  which  solids  are  removed  from  an  eroding  dam.  As  may  be  seen  in 
Table  2,  it  is  common  to  find  expressions  such  as  "dirty  leakage”,  or 
"from  clear  to  chocolate  milk"  as  the  only  available  data  describing  the 


amount  of  solids  being  washed  out  from  the  core  of  a  dam.  Also,  there 
is  very  little  information  on  the  size  of  the  erosion  channel. 

6.  The  authors  recognize  that  there  is  great  difficulty  in  (1)  ob¬ 
taining  measurements  of  the  amount  of  solids  being  eroded  from  a  dam  and 
(2)  obtaining  the  size  of  an  erosion  channel.  However,  the  availability 
of  this  type  of  quantitative  information  is  critical  if  we  are  to  im¬ 
prove  our  present  design  criteria  to  prevent  erosion  from  occurring.  The 
following  paragraphs  summarize  the  current  design  philosophy  on  cracking, 
erosion,  and  the  protection  of  embankment  dam  core  materials. 

7.  Cracking  and  erosion:  Based  upon  the  field  observations  men¬ 

tioned  above,  investigators  have  concluded  that  it  is  possible  thAt  core 
materials  will  crack  regardless  of  design  provisions.  In  general, 
cracks  are  caused  by  differential  settlement  between  the  compressible 
core  and  the  stiffer  shell  material  (load  transfer)  and/or  the  opening 
of  a  crack  due  to  the  fluid  pressure  of  the  reservoir  water  (hydraulic 
fracture) . 

6.  It  is  also  currently  assumed  that  for  cohesionless  materials 
(such  as  filter  materials)  cracks  will  collapse  upon  saturation.  In  co¬ 
hesive  materials,  however,  the  cracks  may  remain  open  due  to  the  un¬ 
drained  strength  of  the  material  (Vaughan,  1976).  Pressurized  water 
from  the  reservoir  may  flow  through  these  cracks  and  cause  concentrated 
leaks  to  appear  downstream.  It  is  this  concentrated  flow  which  has  led 
to  the  erosion  of  core  and  transition  materials  in  a  number  of  dams.  At 
present,  the  "safe"  assumptions  to  make  in  the  design  of  an  embankment 
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dam  are:  "Regardless  of  the  design  considerations  involved,  it  is  pos¬ 

sible  that  the  embankment  dam  material  will  crack"  (Sherard,  1973),  and 
"...  some  erosion  will  occur  from  the  walls  of  a  crack  with  water  flow¬ 
ing  through  it..."  (Vaughan,  1976).  Therefore,  a  filter  is  assumed  to 
be  necessary  to  protect  core  material.  Presently,  filter  design  cri¬ 
teria  are  based  on  the.  grain-size  characteristics  of  core,  filter  and 
shell  materials.  Table  3  presents  a  list  of  design  criteria  as  summa¬ 
rized  by  Hsu  (1981).  Filters,  however,  are  costly  to  install  and  main¬ 
tain  which  makes  their  use  oftentimes  unfeasible,  especially  for  small 
dams  such  as  those  associated  with  small  agricultural  projects. 

Requirements  of  a  Successful  Erosion  Test  Methodology 

9.  Based  upon  field  observations  and  the  present  state  of  knowledge 
of  erosion  in  dams,  investigators  have  concluded  that  the  following  fac¬ 
tors  correlate  with  the  incidence  of  erosion  failures  (Aitchison  and 
Voods,  1965;  Sherard,  1976;  Parkin,  1976): 

•  Rapid  reservoir  filling 

•  The  presence  of  cracks  and  discontinuities 

•  Compaction  on  the  dry  side  of  optimum 

•  Low  ion  concentration  of  the  eroding  fluid 

•  A  high  concentration  of  sodium  ions  in  the  soil  pore  water 

10.  Any  laboratory  test  must  take  into  account  these  factors  as  well 
as  model  1)  the  state  of  stress,  and  2)  the  hydraulic  conditions  under 
which  embankment  dam  erosion  occurs. 


11.  Further,  an  erosion  test  should  be  able  to  model  the  behavior  of 
flow  through  cracked  core  material  into  soil  filters.  Such  tests  are 
important  because  filters  designed  using  the  criteria  presented  in  Table  3  for 
intact  core  material  have  sometimes  failed  to  prevent  erosion  damage  when 
cracks  have  developed.  These  filter  criteria  presented  in  Table  3  were  gener¬ 
ally  developed  from  tests  on  uncracked  base  materials.  Therefore,  the  need  to 
develop  filter  criteria  based  on  laboratory  tests  which  use  soils  with 
open  cracks  has  been  recognized  (Sherard,  1977;  Vaughan,  1978;  Gori  bel* 

lo  and  Lei  Lam,  1979;  Ofegbu,  1981). 

Limitations  of  Existing  Test  Methods 

12.  The  previous  discussion  suggests  that  the  following  two  charac¬ 
teristics  must  be  combined  to  properly  model  the  erosion  behavior  of 
cracked  embankment  core  materials  in  the  laboratory: 

•  The  physical  and  chemical  characteristics  of  the  soil  and  of 
the  eroding  fiuid. 

•  The  insitu  state  of  stress  and  hydraulic  conditions  acting 
on  the  cracked  soil. 

13.  An  indepth  review  of  equipment  and  procedures  available  to  study 
tiie  erosion  potential  of  clays,  prepared  by  Perry  (1975),  provided  in¬ 
formation  that  showed  that  few  if  any  current  testing  procedures  suc¬ 
cessfully  combine  both  of  these  testing  requirements.  In  general,  this 
review  showed  that  there  are  basically  four  classes  of  test  equipment 
and  methodologies  used  to  measure  the  erosion  potential  of  soil  materi¬ 
als.  These  include:  li  the  pinhole  test  and  the  similar,  physical  ero- 
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sion  test  (PET),  2)  the  rotating  cylinder  test,  3)  the  flume  test,  and 
4)  the  hollow  cylinder  test.  These  existing  procedures  do  not  adequate¬ 
ly  model  the  insitu  state  of  stress,  the  effectiveness  of  filters,  or 
the  mechanisms  of  flow  through  a  crack  as  described  in  Table  4. 

Purpose  of  the  Triaxial  Erosion  Test 

14.  To  minimize  the  limitations  discussed  previously,  a  triaxial 
erosion  test  was  developed  by  the  authors  (Sanchez  and  Silver,  1982) 
based  on  the  physical  model  presented  below.  Subsequent  sections  de¬ 
scribe  this  laboratory  triaxial  erosion  test  methodology  and  present  the 
results  of  tests  performed  on  typical  embankment  dam  materials  using 
this  test  method. 


Model  of  the  Physical  Problem 

15.  The  physical  conditions  it.  an  embankment  dam  during  first  fill¬ 
ing  are  presented  diagrammaticallv  in  Figure  la.  If  a  transverse  crack 
exists  in  the  core,  it  will  be  under  a  state  of  stress  as  shown  in  Fig¬ 
ure  lb.  The  total  vertical  stress,  o^.,  will  tend  to  close  the  crack 

while  the  total  horizontal  streaa,  ,  will  tend  to  keep  the  crack  open. 

The  magnitude  of  both  the  vertical  and  horizontal  stress  is  influenced 
by  load  transfer  (arching)  caused  by  differential  settlement  between  the 
shell  and  the  core.  The  compacted  soil  element  is  initially  in  a  par¬ 
tially  saturated  state  under  these  total  stress  conditions. 

16.  As  the  reservoir  is  filled,  a  water  force,  u,  will  act  on  the 
upstream  side  of  the  crack.  The  water  pressure  acting  on  the  downstream 


side  of  the  crack  will  be  zero.  The  resulting  pressure  head,  h,  along 
the  length  of  the  crack,  1  ,  and  the  gradient  ,  h/le,  will  tend  to 

keep  the  crack  open  and  cause  water  to  flow  through  the  crack.  It  is 

this  flow  that  can  erode  material  from  the  surface  of  the  crack. 

\ 

17.  Based  on  the  physical  characteristics  described  above,  the  main 
parameter  required  for  the  study  of  erosion  from  a  cracked  core  material 
is  the  erosion  rate,  5  (defined  as  the  weight  lost  per  unit  area  of 
crack  per  unit  time)  as  a  function  of  the  fluid  shearing  stress,  T,  ap¬ 
plied  to  the  surface  of  the  crack.  Further,  to  better  understand  why 
soils  erode  and  to  be  able  to  quantitatively  determine  the  amount  of  ma¬ 
terial  that  will  be  eroded  from  a  dam,  it  is  important  to  understand  how 
the  erosion  rate,  i,  is  affected  by  changes  in  the  following  variables: 

•  soil  compaction  water  content  and  dry  density 

•  crack  dimensions,  surface  roughness,  and  tortuosity 

•  confining  pressure, 

•  eroding  fluid  gradients,  and 

•  the  chemistry  both  of  the  pore  and  eroding  fluid 

18.  Two  of  these  factors  were  investigated  using  the  test  equip¬ 
ment  and  procedures  described  in  detail  elsewhere  (Sanchez  and  Silver, 
1982)  and  summarized  in  the  following  sections. 
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PART  II:  TRIAXIAL  EROSION  TEST  METHODOLOGY 


Description  of  the  Testing  Equipment 

19.  The  components  of  the  testing  equipment  are  pictured  and  sche¬ 
matically  represented  in  Figures  2  and  3.  The  following  paragraphs 
briefly  describe  the  triaxial  erosion  test  equipment: 

20.  Triaxial  Cell :  The  test  is  performed  in  a  standard  triaxial 
cell  with  art  enlarged  drainage  system  and  modified  top  and  bottom  pla¬ 
tens.  Pictures  of  the  top  and  bottom  platens  are  provided  in  Figure  4. 

21.  Water  Supply  System:  Three  8  liter  (  2  gal .)* capacity  cylindri¬ 
cal  pressure  tanks  were  independently  controlled  to  provide  a  continuous 
high  rate  of  flow  to  the  triaxial  specimen  at  a  desired  constant  pres¬ 
sure  head. 

22.  Effluent  Beakers :  Effluent  containing  eroded  particles  flowing 
out  of  the  specimen  was  collected  and  volumetrically  measured  in  2000  ml 
and  4000  ml  capacity  graduated  beakers.  The  weight  of  the  solids  eroded 
was  determined  by  decanting  the  fluid  from  the  beakers  and  then  drying 
and  weighing  the  solid  residue.  A  screen  equivalent  in  size  to  a  No. 10 
sieve  was  mounted  on  the  bottom  platen  of  the  triaxial  cell  to  help  sup¬ 
port  the  specimen.  The  maximum  particle  size  of  all  specimens  tested 
was  smaller  than  the  No. 10  sieve,  so  therefore  the  bottom  platen  screen 
did  not  interfere  with  the  erosion  of  particles. 

23.  Compaction  Mold:  A  split  compaction  mold  7.1  cm  in  diameter  and 
15.5  cm  in  height  was  used  to  form  specimens.  A  forcing  nozzle,  having 

*  A  table  of  factors  for  converting  U.  S.  customary  unite  of  neasureaent 
to  aetrlc  (SI)  units  la  presented  on  page  5. 
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the  exact  same  shape  as  the  injecting  flow  nozzle  in  the  top  platen 
(Fig. 4),  was  attached  to  the  bottom  plate  of  the  mold.  This  bottom 
plate  and  the  forming  nozzle  were  both  slotted  to  permit  insertion  of  a 
thin  blade  (2.3  cm  wide  by  0.2  cm  thick)  to  form  the  crack.  This  forming  blade 
was  mounted  on  the  head  of  a  specimen  extruding  jack  as  shown  in  Figure  5:  the 
blade  was  jacked  up  slowly  in  increments,  to  avoid  splitting  the  specimen, 
through  the  slotted  bottom  plate  of  the  mold  and  through  the  specimen.  By  this 
process  an  erosion  channel  having  known  dimensions  was  formed  in  the  specimen. 

Test  Procedure 

24.  The  tests  were  performed  using  standard  triaxial  specimens  which 
were  compacted  in  the  special  mold  described  above.  A  triangular  orif¬ 
ice  was  molded  into  one  end  of  the  specimen  by  the  forming  nozzle.  This 
orifice  directed  the  flow  of  water  through  the  specimen  and  prevented 
the  flow  of  water  around  the  circumference  of  the  specimen  or  between 
the  specimen  and  the  confining  membrane.  The  compaction  procedure,  de¬ 
scribed  in  detail  by  Silver  (1976).  was  used  to  minimize  specimen  prepa¬ 
ration  effects  on  soil  behavior;  the  method  consists  of  tamping  moist 
soil  in  layers  (3  layers  for  this  investigation)  with  each  layer  com¬ 
pacted  to  a  prescribed  dry  density.  The  procedure  makes  corrections  for 
the  fact  that  in  compacting  a  material  in  layers,  each  succeeding  laye*- 
densifies  the  material  in  the  layers  below. 

2*.  After  all  the  layers  were  compacted, the  forming  blade  was 
jacked-up  through  the  specimen  to  preform  a  flow  channel  for  the  eroding 
fluid.  The  specimen  was  then  placed  inside  the  triaxial  cell  and  the 
specimen  was  surrounded  by  a  triaxial  membrane  so  that  both  axial  and 
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lateral  confining  pressures  could  be  applied  to  simulate  an  insitu 

stress  condition  (Fig.  1).  Hydraulic  flow  was  modeled  by  allowing  water 
to  flow  from  the  top  to  the  bottom  of  the  specimen.  Distilled  water  was  used 
for  all  erosion  tests  except  for  a  limited  number  of  tests  conducted  to  study 
the  effect  of  eroding  fluid  concentration  (see  paragraphs  70-71).  The  hydrau¬ 
lic  pressure  was  held  constant  at  the  top  of  the  specimen  using  the  large 
volume,  pressurized  reservoirs  described  previously;  the  hydraulic  pressure 
at  the  end  of  the  outflow  line  was  maintained  at  atmospheric  pressure  (Fig.  3). 
This  laboratory  test  therefore  modeled  both  the  stress  and  flow  conditions 
presented  in  Figure  1. 


PART  III:  TEST  PROGRAM 


Objectives  of  this  Study 

26.  The  need  to  understand  how  various  factors  contribute  to  internal 
erosion  or  piping  failure  of  a  cracked  core  in  an  embankment  dam  has 
long  been  recognized.  However,  many  independent  and  related  factors  in¬ 
fluence  erosion  test  results.  Therefore,  the  scope  of  this  report  was 
limited  to  an  evaluation  of  the  influence  of  the  following  important 
factors  within  the  range  of  values  given  in  Table  5: 

•  Compaction  water  content  and  dry  density 

•  Ionic  concentration  of  the  eroding  fluid 

27.  To  minimize  the  number  of  variables,  the  initial  slot  dimensions 
of  2.3  cm  wide  by  0.2  cm  thick,  a  confining  cell  pressure  of  98  kPa 
(14.2  psi),  and  a  hydraulic  head  of  12  kPa  (1.75  psil  equal  to  a  hydrau¬ 
lic  gradient  of  10,  were  all  kept  constant  for  this  investigation.  The 
effect  of  hydraulic  gradient  on  erosion  values  measured  in  the  test  was 
described  In  detail  by  Sanches  and  Silver  (1982).  Varying  the  hydraulic 
gradient  over  the  range  17  to  75  was  found  to  have  little  effect  on  the  erosion 
rate  of  aaterials  of  the  type  studied  in  this  test  program. 

26.  In  addition,  this  investigation  studied  the  validity  of  assump¬ 
tions  made  to  model  the  physical  processes  occuring  in  the  triaxial  ero¬ 
sion  test.  The  investigation  also  determined  the  possibility  of  using 
this  test  method  in  evaluating  the  performance  of  filters  for  the  pro¬ 
tection  of  cracked  fine-grained  soils. 
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Description  of  the  Materials  Tested 


29.  The  study  consisted  of  a  series  of  triaxial  erosion  tests  per* 
formed  on  five  embankment  dam  core  materials,  designated  A,6,C,D,  and  £. 
These  materials  represent  a  wide  range  of  soil  types,  from  a  moderately 
plastic  clay  to  a  silt,  as  may  be  seem  from  their  index  properties  list¬ 
ed  in  Table  6.  Material  E  is  a  slightly  clayey  silt  from  the  core  of 
the  Teton  Dam.  Figures  6  through  10  show  the  grain  size  and  compaction 
characteristics  for  these  materials. 

30.  Recent  research  (Arulanandan  et.al.,  1976)  suggests  that  mechan¬ 
ical  index  properties,  such  as  those  given  in  Table  6,  seem  to  have  lit¬ 
tle  or  no  relation  to  the  erosion  behavior  of  soils.  However,  it  has 
been  extensively  reported  in  both  the  Soil  Science  and  Geotechnical  En¬ 
gineering  literature  that  the  chemistry  of  the  soil  pore  water  may  give 
a  better  indication  of  the  soil  erosion  characteristics.  Table  7  pres¬ 
ents  the  results  of  pore  water  chemistry  tests  for  Materials  A  through 
E.  These  results  can  be  used  to  compare  the  qualitative  erosion  behavior 
predicted  by  identification  charts  given  in  the  literature  (e.g.  Sherard 
et.al.,  1976)  to  the  quantitative  values  obtained  in  this  investigation. 
Appendix  C  describes  the  procedure  for  obtaining  soil  pore  water  ex¬ 
tracts  . 
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PART  IV:  METHOD  OF  DATA  COLLECTION 


31.  Data  obtained  for  each  test  included  1)  a  description  of  the 
physical  characteristics  of  the  erosion  test  specimen  both  before  and 
after  the  test,  2)  the  hydraulic  parameters  of  the  test  system,  3)  the 
quantity  of  fluid  flowing  through  the  specimen  as  a  function  of  time, 
and  4)  the  amount  of  solid  material  eroded  from  the  specimen  as  a  func¬ 
tion  of  time.  Methods  used  to  obtain  and  calculate  these  values  are  de¬ 
scribed  in  more  detail  in  the  following  paragraphs. 

Specimen  Characteristics 

32.  Atterberg  Limit  tests,  hydrometer  tests  and  compaction  tests  were 
performed  on  all  materials  tested.  Initial  specimen  density  was  deter¬ 
mined  using  measurements  taken  with  a  Pi-tape  and  a  scale,  and  by  weigh¬ 
ing  the  specimen.  Specimen  height  changes  were  monitored  during  each 
erosion  test  using  a  dial  gage  attached  to  the  loading  rod  of  the  triax- 
ial  cell.  After  the  test,  the  shape  of  the  eroded  crack  was  qualita¬ 
tively  described.  The  dimensions  of  the  specimen  were  again  measured 
and  the  specimen  was  oven  dried  and  weighed  to  determine  its  final 
weight. 

Erosion  Test  Parameters 

33.  Test  parameters  measured  and  monitored  during  the  erosion  test 
included:  1)  fluid  pressure  at  the  top  of  the  specimen;  2)  volume  of 
eroded  fluid  and  solids  as  functions  of  time;  and  3)  height  change  of 
the  specimen.  Both  cell  pressure  and  head  water  pressure  were  measured 
to  an  accuracy  of  5  kPa  (0.7  psi).  The  volume  of  the  eroded  fluid  was 
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recorded  to  an  accuracy  of  50  ml.  Time  was  measured  to  the  nearest  0.01 
min  and  dry  weight  of  the  eroded  solids  was  recorded  to  the  nearest  0.01 
gram. 

34.  Prior  to  the  test,  the  confining  cell  pressure  and  water  supply 
head  pressure  were  read  with  a  calibrated  gage  or  pressure  transducer. 
During  the  test,  flow  into  the  effluent  receiving  beakers  was  continu¬ 
ously  timed  at  flow  intervals  of  between  400  and  500  ml.  The  time  was 
recorded  when  a  beaker  was  filled  to  capacity  with  fluid  and  eroded  ma¬ 
terial,  and  the  next  beaker  began  to  fill. 

35.  The  test  parameters  described  above  were  used  to  determine  the 
relationship  between  material  properties,  hydraulic  parameters,  and  the 
time  rate  of  material  erosion.  Typical  data  showing  these  relationships 
is  presented  in  subsequent  sections. 


PART  V:  METHOD  OF  DATA  ANALYSIS 


36.  Data  from  each  erosion  test  was  analyzed  to  obtain  values  of  the  ero¬ 
sion  rate,  c  (defined  as  the  weight  lost  per  unit  cross-sectional  area  of  the 
crack  per  unit  time)  as  a  function  of  the  eroding  fluid  shear  stress,  t  .  To 
obtain  these  parameters,  calculations  were  first  made  to  determine  the  varia¬ 
tions  of  erosion  weight  loss  versus  time  (g/min)  and  the  quantity  of  flow 

3 

passing  through  the  crack  versus  time  (cm  /min). 

37.  Computation  of  Hydraulic  Gradient :  The  hydraulic  pressure  at 

the  end  of  the  outflow  line  was  maintained  at  atmospheric  pressure. 
Therefore,  the  hydraulic  gradient  was: 

h  ‘  hb  (1) 

i  =  —  ■ — ■— 

1 

e 


where: 

i  *  hydraulic  gradient; 

h  *  head  of  water  applied  to  the  top  of  the  specimen  in  me¬ 
ters; 

h.  *  head  of  water  applied  at  the  end  of  the  outflow  line 

D 

(.equal  to  zero  in  this  study)  ; 

1^  *  eroded  crack  length  in  meters  (obtained  by  subtracting 
the  nozzle  height  from  the  final  height  of  the  specimen). 
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38.  Initial  Crack  Cross-Section:  It  was  assumed  that  1)  the  initial 


cross-sectional  area  of  the  crack,  ,  and  2)  the  cross-sectional 

area  of  the  crack  after  confining  pressure  was  applied  to  the  specimen, 
but  before  initiating  the  flow,  (0),  were  equivalent  and  equal  to 

the  cross-sectional  area  of  the  blade  forming  the  crack  (Fig.  lb): 


A  =  A  (0)  =  d  •  w 
x  x 


(2) 


where : 


A^  =  initial  cross-sectional  area  of  the  crack,  equal  to  the 
cross-sectional  area  of  the  forming  blade; 

A^  (0)  =  crack  cross-sectional  area  at  time  equal  to  zero, 

after  application  of  confining  pressure; 

d  =  thickness  of  the  forming  blade  and  of  the  crack; 

w  =  width  of  the  forming  blade  and  of  the  crack. 

This  assumption  implies  that  there  was  little  if  any  closure  of  the 
crack  as  a  result  of  applying  the  confining  pressure. 

39.  Variation  of  Cross-Sectional  Area  with  Time:  To  compute  erosion 
rates  and  fluid  shear  stresses,  numerical  values  for  the  surface  area 
and  cross-sectiona.  area  of  the  crack  were  required  as  a  function  of 
time.  Assumptions  were  required  to  compute  these  areas  as  no  direct 
measurement  could  be  made  during  the  test.  One  of  these  assumptions  was 
made  based  on  the  observation  that  srter  tne  completion  of  a  test,  the 


erosion  surface  did  not  vary  significantly  along  the  length  of  the 
crack.  For  this  condition,  the  variation  of  the  cross-sectional  area  of 
the  crack,  A^  (t) ,  with  time  can  be  shown  to  be  (Sanchez  and  Silver, 

1982): 


A  (t) 
x 


P  •  t 

-  +  A  (0) 

x 


(3) 


where: 

A  (t)  =  cross -sectional  area  of  crack  as  a  function  of 

x 

time  in  cm*  : 

P  =  slope  of  cumulative  dry  weight  eroded  versus  time  plot  in 
g/min; 

t  =  time  in  minutes; 
p  *  dry  density  of  specimen  in  g/cm*  ; 

ip  *  eroded  length  of  the  specimen  at  the  end  of  the  test  in 
cm;  and 

A  (0)  *  cross-sectional  area  at  t  *  0  in  cm1, 
x 

40.  The  eroded  length  of  the  specimen,  1  (equal  to  the  initial 

© 

height  of  the  specimen  minus  the  nozzle  length  and  the  final  vertical 
dial  gauge  reading;  was  regarded  as  a  constant  since  it  varied  less  than 
about  0.2  cm  during  any  test. 


41.  Variation  of  Erosional  Surface  Area:  After  the  cross-sectional 


area  was  computed,  the  wetted  perimeter  of  the  erosion  channel,  and 
hence  the  erosional  surface  area  was  calculated.  Observations  made  at 
the  end  of  a  test  indicated  that  in  almost  all  cases  the  final  eroded 
shape  was  roughly  elliptical.  However,  for  this  report  it  was  assumed 
that  the  crack  enlarged  in  a  rectangular  fashion,  as  shown  in  Figure  1. 
The  rectangle  was  assumed  to  have  a  constant  width  equal  to  the  width  of 
the  slot-forming  blade;  only  the  thickness  was  assumed  to  vary  with  time 
to  account  for  the  soil  volume  being  eroded.  Using  this  assumption,  the 
erosional  surface  area  as  a  function  of  time,  A^  (t),  can  be  shown  to 

be  equal  to: 


A  (t)  =  2  •  [  A  (t)/w  +  w]  •  1 

5  X  © 


(4) 


where: 


A  ft)  -  erosional  surface  area  of  crack  as  a  function  of 
s 

time  in  cm*  ; 

w  s  width  of  the  forming  blade  and  of  the  crack  m  cm;  and 

1^  2  eroded  length  of  the  specimen  at  the  end  of  test  in  cm. 

42.  In  Appendix  D,  other  assumptions  were  made  for  the  time  rate  of 
change  of  the  shape  of  the  erosion  channel.  Results  of  analyses  assum¬ 
ing  a  constant  width  elliptical  erosion  channel  and  a  variable  width  el¬ 
liptical  erosion  channel  were  compared  to  the  constant  width  rectangular 
shape  assumption  presented  i;:  Equatic::  -.  These  analyses  indicated  that 
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using  the  simple  rectangular  assumption  resulted  in  erosion  channel  size 
parameters  within  ±  15"o  of  the  possibly  more  accurate  elliptical  erosion 
channel  size  parameters.  This  small  difference  was  felt  to  be  insignif¬ 
icant  for  the  purposes  of  this  study. 

43.  Eroding  Fluid  Shear  Stress :  The  flow  through  the  crack  was  ana¬ 
lyzed  as  a  closed  channel  pipe  flow.  Reynold's  numbers  were  calculated 
and  in  all  cases  they  were  less  than  about  2000  indicating  laminar  flow. 
Assuming  that  laminar  flow  did  indeed  occur,  the  fluid  shear  stress,  x, 
was  derived  as  (John  and  Haberman.  1971): 


x 


*  V  y 
Dh  5 


(5) 


where: 

x  =  fluid  shear  stress  in  N/m*  ; 

V  ■  average  fluid  velocity  [  equal  to  Q/Ax  (t))  in  cm/min  ; 

V  *  viscosity  of  water  in  N  •  sec/m*  ; 

*  hydraulic  diameter  of  the  crack  in  cm  ;  equal  to 


£  =  conversion  factor  equal  to  .60  sec/min  ;  and 


«  flow  parameter  (equal  to  12  for  this  study). 
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44.  Appendix  E  presents  detailed  results  of  a  series  of  hydraulic  ex¬ 
periments  performed  to  check  the  validity  of  the  assumption  that  flow 
through  the  erosion  channel  was  laminar.  These  experiments  showed  that 
some  flow  separation  could  be  introduced  into  the  flow  at  the  abrupt 
transition  between  the  narrow  injecting  flow  nozzle  and  the  rapidly  en¬ 
larging  flow  channel. 

45.  An  analysis  of  this  flow  separation  indicated  that  the  equations 
for  laminar  shear  stress  were  still  applicable  if  the  value  of  the  shear 
stress  was  increased  (Olson,  1980).  Thus,  the  value  of  the  flow  parame¬ 
ter,  i>,  in  Equation  5  was  increased  by  30°.  over  its  average  laminar  val¬ 
ue  of  9.3  to  an  adjusted  value  of  12  to  account  for  this  possible  flow 
separation.  Appendix  E  presents  more  details  of  this  analysis. 

46.  Erosion  Rate:  The  erosion  rate,  t,  was  computed  as: 


P  (7) 

l  *  - 

As  (t) 

where: 

c  *  erosion  rate  g/(cra2  •  rain)  ; 

P  *  slope  of  cumulative  weight  eroded  versus  time  curve  (g/ 
min)  : 

(t)  *  erosional  surface  area  of  crack  in  cm1. 
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47.  Solids  Eroded  as  a  Function  of  Volume  of  Flow:  The  literature 


suggests  that  estimates  of  the  amount  of  erosion  occurring  through  the 
core  of  a  dam  may  be  obtained  from  field  measurements  of  the  weight  of 
solids  and  the  volume  of  water  flowing  from  the  dam  (Bertram,  1967; 
Kjaernsli  and  Torblaa,  1968).  For  comparison  with  these  published  val¬ 
ues,  an  erosion  factor,  S,  was  defined  and  calculated  as  the  dry  weight 
of  solids  eroded  divided  by  the  volume  of  water  collected. 

48.  Percent  of  Weight  Erosion  per  Unit  Time:  Another  useful  value 
that  may  be  used  to  compare  the  results  of  this  investigation  to  the  re¬ 
sults  obtained  by  others  is  the  percent  erosion  per  unit  time,  a.  This 
value  is  defined  as  follows: 


P  •  100*,.  (8) 

a  =  —  ■■ .  . 

V. 

l 


where: 


a  =  the  percent  erosion  per  unit  time  in  percent/min  ; 

P  =  slope  of  cumulative  dry  weight  eroded  versus  time  plot  in 
g/min;  and 

-  Dry  weight  of  the  specimen  at  the  beginning  of  a  test 
in  grams. 

It  should  be  noted  that  a  similar  erosion  quantity  is  often  given  in  the 
literature.  That  quantity  is  simply  called  "Percent  Erosion"  and  is  de¬ 
fined  as  the  cumulative  weight  of  material  eroded  divided  by  the  initial 


.  Lewis  and  Schmidt,  1977).  It  may  be 


weight  of  the  specimen,  VT  (e.g 

seen  that  "Percent  Erosion"  does  not  include  the  effect  of 
tion.  It  therefore  has  limited  use  in  quantifying  erosion, 
specimens  are  tested  for  the  same  amount  of  time. 


27 


-hi 


test  dura- 
unless  all 


I 


PART  VI:  PRESENTATION  OF  DATA 


49.  The  results  of  typical  triaxial  erosion  tests  for  specimens  of 
different  embankment  core  material  types  are  shown  in  Figures  11  through 
15.  These  figures  correlate  specimen  characteristics  and  experimental 
values  measured  in  the  erosion  test.  These  experimental  values  include 
the  cumulative  volume  of  flow  and  the  cumulative  weight  of  solids  eroded 
from  the  walls  of  the  crack,  both  as  functions  of  time.  Plotted  in  each 
figure  are  typical  results  for  a  silt  (high  erodibility)  and  for  a  clay 
(low  erodibility)  to  indicate  the  range  of  values  measured  in  the  test¬ 
ing  program. 

50.  For  example,  Figure  11  gives  the  form  of  test  results  showing 
the  volume  of  fluid  flowing  through  a  crack  as  a  function  of  time.  This 
straight  line  relationship,  with  a  small  range  of  variation  in  all 
tests,  indicates  that  the  flow’  rate  is  constant  under  any  applied  con¬ 
stant  gradient.  The  slope  of  this  line  which  is  the  flow  rate,  Q 
(cm’/min),  has  been  computed  by  linear  regression  for  e.ach  test  and  is 
given  in  Table  8. 

51.  Figure  12  shows  typical  results  giving  the  cumulative  dry  weight 
of  solids  eroded  from  the  walls  of  the  crack  as  a  function  of  time. 
This  relationship  in  all  tests  was  also  found  to  be  linear.  The  slope, 
denoted  by  P  (g/rain),  for  each  test  was  computed  and  is  also  shown  in 
Table  8.  Values  for  P  were  high  for  the  silt  core  materials  and  lower 
for  the  clay  core  materials. 
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52.  Typical  plots  of  the  cumulative  dry  weight  of  solids  eroded  as  a 
function  of  volume  of  flow  are  shown  in  Figure  13.  The  slopes  of  these 
lines,  the  erosion  factor,  S  (ppm),  presents  a  useful  quantity  that  has 
been  cited  in  dam  literature  to  provide  a  quantitative  value  for  the 
amount  of  erosion  measured  in  the  field.  It  is  thus  useful  to  relate 
field  measurements,  in  terms  of  embankment  material  characteristics  and 
estimates  of  the  crack  dimensions,  to  results  of  laboratory  tests  to 
provide  criteria  for  the  evaluation  of  cracking  and  erosion  in  existing 
and  future  dams.  These  test  results  found  high  values  of  S  for  the  silt 
core  materials  and  lower  values  of  S  for  clay  core  materials. 

53.  Figure  14  presents  typical  plots  of  cumulative  weight  of  solids 
eroded  per  unit  area  of  crack  surface  as  a  function  of  time.  This  data 
may  be  used  with  Equation  7  to  calculate  the  erosion  rate,  I.  Higher 
erosion  values  were  measured  for  silt  core  materials  than  for  clay  core 
materials. 

54.  Figure  15  plots  typical  curves  of  erosion  rate,  f ,  as  a  function 
of  the  fluid  shear  stress,  t .  Typical  computations  for  these  quantities 
are  shown  in  Table  9.  Note  that  both  the  fluid  shear  stress  and  the 
erosion  rate  were  highest  at  the  beginning  of  a  test  and  decreased  with 
time  which  on  these  plots  increases  towards  the  left  of  the  figure. 

55.  Detailed  test  results  of  the  form  presented  above  for  all  speci¬ 
mens  tested  are  given  in  Appendix  A.  Note  that  the  plots  presented  in 

Appendix  A  that  Involve  measurements  of  the  volume  of  flow  are  drawn  to 
two  different  scales  because  flow  quantities  before  failure  (when  the  speci¬ 
men  caved  in  at  a  severely  eroded  point  on  the  sample  surface)  were  about  ten 

times  smaller  for  silt  specimens  than  for  clay  specimens. 
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56,  Table  8  summarizes  the  results  in  Appendix  A  by  presenting  spec¬ 
imen  measurements  and  test  values  for  the  flow  rate,  Q,  the  rate  of 
weight  erosion,  P,  and  the  ratio  of  solids  eroded  to  the  volume  of  flow, 
S,  at  any  time  t.  It  may  be  noted  that  these  values  of  Q,  P,  and  S  were 
obtained  without  making  any  assumption  about  the  shape  of  the  erosion 
channel.  Therefore,  these  values  provide  comparative  data  with  which  to 
relate  results  obtained  from  this  test  method  with  results  obtained  from 
other  laboratory  and  field  test  methods  that  do  not  include  measurements 
of  the  erosional  surface  area. 

57.  An  interesting  hypothesis  may  be  made  based  upon  the  results  giv¬ 
en  in  Appendix  A  by  noting  that  the  erosion  rate  asymptotically  ap¬ 
proached  a  constant  value  as  the  fluid  shear  stress  increased.  There¬ 
fore,  the  erosion  rate  at  which  the  slope  of  the  erosion  rate  versus 
fluid  shear  stress  curve  was  a  minimum,  was  defined  as  "the  maximum  ero¬ 
sion  rate"  for  each  specimen.  It  should  be  noted  that  the  maximum  ero¬ 
sion  rate  will  vary  from  one  specimen  to  another  depending  on  the  com¬ 
paction  and  eroding  fluid  used  to  test  the  material.  Typical  values  of 

the  maximum  erosion  rate,  t  ,  are  shown  in  Figure  15  and  given  in 

max 

Table  8.  It  should  also  be  noted  that  this  maximum  erosion  rate  value 
is  a  simple  constant  which  takes  into  account  the  erosional  surface 
area.  Therefore,  the  maximum  erosion  rate  provides  a  better  method  of 
quantifying  erosion  chan  other  parameters  which  do  not  include  the  ero¬ 
sional  surface  area  such  as  percent  erosion,  percent  erosion  per  unit 
time  (a),  and  tiie  ratio  of  solids  to  volume  of  fluid  factor  (S).  Appen¬ 
dix  F  presents  a  possible  method  in  which  the  maximum  erosion  rate  pa¬ 
rameter  may  be  used  for  design. 
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PART  VII:  TEST  RESULTS  ON  VARIABLES  STUDIED 


58.  The  effect  of  material  and  fluid  parameters  (including  molding 
water  content,  dry  density,  and  eroding  fluid  ionic  concentration)  on 
the  erosion  potential  of  a  wide  range  of  embankment  core  material  types 
was  investigated.  Detailed  test  results  are  reported  for  Material  B 
which  was  a  slightly  silty  clay,  and  for  materials  D  and  E  which  were 
slightly  clayey  silts. 

59.  Two  more  materials,  A  and  C,  were  also  tested  in  this  experimen¬ 
tal  program.  The  results  of  tests  performed  on  Materials  A  and  C  are 
presented  in  Table  10.  Material  A  was  a  clay,  while  Material  C  was  a 
sandy  clay.  When  tested,  Material  A  exhibited  extremely  high  resistance 
to  erosion.  For  example,  after  5  hours  in  the  triaxial  erosion  test  ap¬ 
paratus,  less  than  7  grams  of  solids  were  eroded  from  Material  A  speci¬ 
mens.  On  the  other  hand,  for  materials  B,  D,  and  E,  between  35  and  133 
grams  of  solids  were  eroded  from  the  specimens  within  30  minutes  after 
the  beginning  of  a  test.  This  large  difference  in  the  erosion  charac¬ 
teristics  of  Material  A  as  compared  to  the  other  materials  tested  may  be 
also  observed  by  comparing  the  maximum  erosion  rate  values  measured  for 
each  material.  The  maximum  erosion  rate  for  Material  A  was  only  0.0004 
g/min-cm*  ,  whereas  the  maximum  erosion  rate  values  ranged  from  a  mini¬ 
mum  of  0.01  g/min-cm*  for  Material  B  to  a  maximum  value  of  about  0.35 
g/min-cm*  for  Material  D  or  E.  Thus,  further  tests  with  Material  A  were 
not  performed. 

*>0.  Tests  were  also  performed  on  a  sandy  clay.  Material  C,  but  the 
flow  rate  during  a  typical  test  decreased  rapidly  from  the  usual  average 
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of  about  500  cm1 /min  to  a  low  value  of  only  100  cm1 /min.  Almost  no 
weight  erosion  took  place  at  this  very  low  flow  rate.  It  was  observed 
that  the  high  percentage  of  sand  particles  (32%  coarser  than  the  No. 200 
sieve)  in  this  material  was  causing  the  crack  to  clog  up.  Therefore, 
additional  tests  were  performed  on  Material  C  to  study  the  effect  of  re¬ 
moving  some  sand  from  the  material.  It  was  found  that  the  flow  rate  and 
thus  erosion  increased  when  material  greater  than  the  No. 50  «ieve  (23% 
retained)  was  removed. 

61.  Further  tests  are  needed  to  determine  whether  the  sand  particles 
in  Material  C  actually  caused  the  crack  to  seal  internally,  or  whether 
the  sand  particles  were  being  lodged  on  dead  areas  of  the  bottom  platen 
of  the  apparatus,  as  described  in  Appendix  E. 

62.  The  use  of  an  upstream  sandy  filter  has  been  recommended  by  dam 
consultants  for  many  years  (Wilson  and  Marsal,  1979).  Test  results  on 
Material  C  suggests  that  such  filters  might  have  been  effective  in  pre¬ 
venting  total  failure  in  danr>  which  had  cracking  problems.  The  behavior 
of  Material  C  suggests  that  if  enough  sand  particles  are  washed  into  a 
crack,  the  crack  will  clog  resulting  in  reduced  seepage  and  erosion. 

63.  The  erosion  rate  versus  shear  stress  curves  presented  in  Appen¬ 
dix  A  have  been  summarized  in  Appendix  B  for  each  material  tested  tc 
show  the  effects  of  molding  water  content,  dry  density,  and  eroding  flu¬ 
id  ionic  concentration  on  erosion.  A  listing  of  the  specimens  tested 
and  additional  computations  performed  to  observe  these  effects  have  also 
been  provided  in  Appendix  B. 


64.  Molding  Water  Content:  From  the  data  in  Appendix  B  the  effect 


of  specimen  molding  water  content  on  the  erosion  characteristics  of  typ¬ 
ical  embankment  core  materials  was  investigated. 

65.  Figure  16  shows  the  effect  of  molding  water  content  on  the  maxi¬ 
mum  erosion  rate,  e  ,  for  Materials  B,  D,  and  E.  Figure  17  shows 

*  max  *  *  ’  ° 

the  effect  of  molding  water  content  on  another  erosion  parameter,  per¬ 
cent  erosion  per  unit  time  (  a  ).  This  quantity  provides  comparative 
data  between  the  triaxial  erosion  test  and  other  laboratory  test  methods 
that  do  not  take  into  account  the  erosional  surface  area. 

66.  A  comparison  of  the  data  presented  in  Figures  16  and  17  shows 
that  the  molding  water  content  has  a  large  influence  on  the  erosion  of 
silt  materials  (Materials  D  and  E).  This  effect,  however,  is  less  pro¬ 
nounced  for  clay  sMterials  (Material  B).  These  results  suggest  that  close 
control  of  the  compaction  water  content  during  construction  (if  feasible  from 

a  practical  viewpoint)  may  minimize  the  erosion  potential  of  silt  but  not  clay. 

67.  Dry  Density:  The  effect  of  specimen  dry  density  on  the  erodi- 

bility  of  typical  embankment  core  materials  was  also  investigated  using 
the  data  in  Appendix  B. 

66.  Figure  18  shows  the  effect  of  specimen  dry  density  on  the  maxi¬ 
mum  erosion  rate,  c  ,  for  Meterials  B,  D,  and  E  compacted  at  the 

RAX 

optimum  moistute  content.  Another  indication  of  the  effect  of  specimen 
dry  density  on  erosion  for  these  materials  is  presented  on  Figure  14 
which  plots  percent  erosion  per  unit  time,  a. 
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69.  Both  of  these  figures  show  that  there  is  little  effect  of  dry 
density  on  erosion  for  silt  specimens  compacted  at  the  optimum  molding 
water  content.  However,  the  clay  material  (Material  B)  was  more  eroda- 
ble  as  the  density  decreased. 

70.  Eroding  Fluid  Concentration:  The  effect  of  the  ionic  concentra¬ 
tion  of  the  eroding  fluid  on  the  erodibility  of  a  clay,  Material  b,  and 
of  a  silt,  Material  E,  was  also  investigated.  Two  ion  concentrations 
were  used  (3  meq/1  and  6  meq/1).  These  concentrations  values  are  within 
the  range  of  typical  reservoir  waters.  The  3  meq/1  solution  used  in  the 
study  was  tap  water  with  2/3  calcium  ions,  1/3  magnesium  ions  and  a 
trace  of  sodium.  The  6  meq/1  concentration  was  prepared  on  the  basis  of 
one  third  equal  parts  of  sodium,  magnesium  and  calcium  ions. 

71.  Figure  20  shows  that  the  maximum  erosion  rate,  t  ,  decreased 

®  max 

slightly  (as  expected)  when  the  eroding  fluid  ionic  concentration  was 
increased  from  3  meq/1  to  6  meq/1.  This  effect  was  slightly  more  pro¬ 
nounced  for  silt  than  for  clay. 

72.  Use  of  Standard  Concrete  Sand  as  Filter  Material :  The  possibil¬ 
ities  of  applying  the  triaxial  erosion  test  methodology  to  test  the  ef¬ 
fectiveness  of  filter  materials  for  the  protection  of  cracked  fine¬ 
grained  soils  was  investigated. 

73.  Figure  21  shows  the  grain  size  limits  for  a  standard  concrete 
sand  material  (ASTM  C - 33 ) .  This  figure  also  shows  the  grain  size  dis¬ 
tribution  of  the  standard  concrete  sand  used  as  the  filter  for  this  in¬ 
vestigation,  and  tne  grain  size  distribution  of  a  typical  base  soil  (Ma¬ 
terial  Ei. 
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,4.  The  filter  materiel  mas  compacted  in  three  layers  to  a  dry  den¬ 
sity  of  1.763  Mg/m>  (110  lb/ft' 1  at  »  mater  content  of  124  and  then  fro- 
gase  specimens  mere  prepared  in  the  same  manner  as  for  the  other 
erosion  test  specimens.  The  base  specimen  mas  then  placed  on  top  of  the 
filter  specimen  in  the  triaxi.l  erosion  cell  and  the  filter  mas  .Homed 
to  them  for  about  two  hoars  prior  to  testing.  The  direction  of  flow  in 
the  test  mas  domnm.rd  through  the  cracked  base  material  and  into  the  in- 

tact  sand  filter. 

The  results  of  these  tests,  summarised  in  Table  11.  shomed  that 
the  standard  concrete  sand  mas  effective  in  eliminating  the  loss  of  base 
material  of  even  the  highly  erosive  Material  E  (Teton  Dam  material,  at 
hydraulic  heads  of  up  to  40  meters  which  mas  equivalent  to  a  gradient  of 
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PART  VIII:  DISCUSSION  OF  THE  TEST  RESULTS 


76.  Results  of  this  study  were  compared  to  results  obtained  by  other 
researchers.  The  following  paragraphs  describe  these  comparisons  and 
discuss  the  test  results  in  more  detail. 

77.  Effect  of  Molding  Water  Content :  Figures  16  and  17,  summarizing 
the  effect  of  molding  water  content  on  erosion,  show  an  initial  de¬ 
crease,  an  inflection  point,  and  then  an  increase  in  both  the  maximum 

erosion  rate,  e  ,  and  percent  erosion  per  unit  time,  a,  as  the  wa* 
max 

ter  content  increases.  Similar  behavior  was  observed  and  reported  by 
Lewis  and  Schmidt  (1977)  for  tests  performed  on  an  Oklahoma  silty  clay 
using  the  pinhole  test.  The  plastic  and  liquid  limits  for  their  materi¬ 
al  were  given  as  15*.  and  31%,  respectively.  Their  results  are  presented 
in  Figure  22.  The  results  of  the.ir  investigation  and  the  results  ob¬ 
tained  in  this  report  for  a  similar  material  (Material  B)  have  also  been 
replotted  in  terms  of  percent  erosion  per  unit  time,  a,  as  shown  in  Fig¬ 
ure  23. 

76.  Another  similar  study,  conducted  by  Grissinger  (1962),  used  the 
flume  test  to  study  the  effect  of  compaction  water  content  on  the  ero¬ 
sion  potential  of  a  Grenada  silt  loam  (Perry,  1975).  The  plastic  and 
liquid  limits  for  this  material  were  given  as  20  and  31  percent,  respec¬ 
tively.  The  optimum  water  content  was  not  provided,  but  an  estimated 
value  of  21?.  was  calculated  from  the  Atterberg  Limits  (Bowles.  1970). 
The  results  of  Grissinger’s  investigation  have  been  replotted  in  terms 
of  the  maximum  erosion  rate  and  are  shown  in  Figure  24.  This  figure 
shows  that  the  erosion  behavior  of  the  Grenada  silt  loam  was  similar  to 
the  behavior  measured  for  Materials  B,  D,  and  E. 
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79.  Lewis  and  Schmidt  noted  that  the  initial  decrease  in  erosion  with  in¬ 


creasing  water  content  was  contrary  to  expectations  from  a  consideration  of  the 
theory  of  structure-stability  of  clay  particles.  This  theory  suggests  that 
there  are  two  principal  forces  which  act  between  adjacent  particles:  attrac¬ 
tive  (or  Van  der  Wall)  forces  between  the  atoms  of  different  particles,  and 
repulsive  (or  Osmotic)  forces  between  the  electrical  double  layer  sur¬ 
rounding  each  particle.  Higher  compaction  water  content  results  in  a 
larger  distance  between  the  clay  particles.  Therefore,  since  the  attrac¬ 
tive  Van  der  Wall  forces  are  inversely  proportional  to  some  high  power 
of  the  interparticle  distance,  the  attractive  forces  decrease  with  in¬ 
creasing  water  content  due  to  the  increased  particle  distance.  Thus, 
the  erosion  should  have  increased  as  the  water  content  increased. 

80.  In  addition,  higher  compaction  water  content  should  yield  a  lower 
ionic  concentration.  This  lower  ionic  concentration  has  been  reported 
to  produce  a  stronger  repulsion,  or  an  increase  in  osmotic  forces,  be¬ 
tween  particles  (Resendiz,  1977).  Therefore,  this  lower  ionic  concen¬ 
tration  should  have  also  resulted  in  an  increase  in  erosion.  This  dis¬ 
cussion  may  explain  the  behavior  shown  in  Figures  16,  17,  22,  23,  and  24 
for  water  contents  above  optimum,  but  it  does  not  explain  the  decrease 
in  erosion  observed  for  water  contents  below  optimum. 

81.  The  initial  decrease  in  erosion  for  water  contents  below  optimum 
was  attributed  by  Lewis  and  Schmidt  to  slaking.  Slaking  is  the  break¬ 
down  of  soil  upon  immersion  in  a  fluid.  One  of  the  mechanisms  that  con¬ 
tributes  to  this  soil  breakdown  is  the  action  of  capillary  absorption  of 
the  fluid  and  its  subsequent  compression  of  the  air  within  the  voids  of 
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the  soil.  This  compression  may  become  higher  than  the  forces  keeping 
the  soil  particles  together,  causing  slaking.  A  reduction  in  slaking 
and  therefore  in  erosion  should  occur  as  the  amount  of  entrapped  air 
available  for  compression  decreases  as  the  water  content  increases. 

82.  The  inflection  point,  shown  in  Figures  16,  17,  and  22,  is  of 
particular  interest  since  it  corresponded  to  the  point  of  minimal  ero¬ 
sion.  For  the  tests  performed  by  Lewis  and  Schmidt  on  the  Oklahoma  sil¬ 
ty  clay,  this  inflection  point  coincided  with  the  water  content  at  both 
the  plastic  limit  (15%)  and  the  standard  Proctor  optimum  moisture  con¬ 
tent  (14%).  Therefore,  the  change  from  decreasing  erosion  to  increasing 
erosion  was  attributed  to  the  amount  of  water  within  the  voids  at  or 
above  the  plastic  limit. 

83.  It  may  be  seen  that  Grissinger's  results  are  of  the  same  form, 
and  that  the  inflection  point  also  coincides  with  the  water  content  at 
both  the  plastic  limit  (20%)  and  at  the  standard  Proctor  optimum  mois¬ 
ture  content  (21%).  However,  Lewis  and  Schmidt's  results  for  specimens 
compacted  by  the  modified  Proctor  method,  as  shown  in  curve  B  of  Figure 
23,  showed  the  inflection  point  being  somewhat  shifted  towards  the  modi¬ 
fied  Proctor  optimum  water  content  of  10%. 

84.  For  the  materials  tested  for  this  report,  however,  the  inflection 
point  did  not  coincide  with  the  plastic  limit;  the  inflection  point  was 
1  or  2  percentage  points  above  the  optimum  water  content  of  each  materi¬ 
al  while  the  plastic  limit  was  about  6  percentage  points  above  the  opti¬ 
mum  water  content  for  all  three  materials.  Thus,  it  appears  that  the 
optimum  compaction  moisture  content  is  a  better  indicator  oi  the  mecha¬ 
nism  affecting  erosion  than  the  plastic  limit. 
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85.  Effect  of' Molding  Water  Content  on  Teton  Dam  Material : 
interesting  to  note  that  the  core  of  the  Teton  dam  embankment  was  com¬ 
pacted  at  an  average  water  content  of  18.5?0,  or  1%  dry  of  the  core's  av¬ 
erage  optimum  water  content  of  19.6%  (U.S  Department  of  the  Interior, 

1977).  It  may  be  seen  in  Figure  16,  that  compacting  Teton  Dam  specimens 
(Material  E)  at  a  molding  water  content  1%  dry  of  optimum,  gave  an  ero¬ 
sion  rate  of  0.3  g/min-cm2  which  was  twice  as  high  as  the  minimum  meas¬ 

ured  erosion  rate  obtained  at  +1.5  %  wet  of  optimum  of  0.16  g/min-cm2. 

Both  of  these  erosion  rates  are  high.  Erosion  tests  which  provide  quantitative 

data  could  be  useful  to  the  dam  designer  faced  with  making  a  choice  between  two 

nondispersive  borrow  materials  for  the  core  of  a  dam  or  a  choice  of  compaction 
water  contents  on  a  given  core  material. 

86.  Effect  of  Density  on  Erosion:  Figures  18  and  19  summarize  the 

effect  of  variation  in  specimen  density  on  erosion.  Limited  data  shows 
that  there  v.as  little  effect  of  specimen  density  on  erosion  from  95?»  to 
99%  of  standard  Proctor  maximum  density  for  the  silt  materials  D  and  E. 
However,  there  was  some  decrease  in  erosion  when  the  density  of  clay 
(Material  B)  was  increased  from  90%  to  98%  of  standard  Proctor  density. 

87.  Figure  25  further  evaluates  the  effect  of  density  on  erosion  by- 
plotting  percent  erosion  per  unit  time  versus  percent  standard  Proctor 
density  for  the  clay  material  B  and  an  Oklahoma  silty  clay  studied  by 
Lewis  and  Schmidt  (1977).  Data  from  Lewis  and  Schmidt  were  replotted 
using  only  specimens  which  were  at  about  the  same  moisture  content  of 
either  13%  to  15%  or  7%  to  8%.  The  higher  moisture  content  corresponded 
with  both  the  plastic  limit  and  the  standard  Proctor  optimum  moisture 
content  of  15%. 
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8f .  No  effect  of  density  on  erosion  was  noted  for  specimens  at  a  wa¬ 
ter  content  of  13%  to  15%  which  was  at  about  the  standard  optimum  and 
the  pi.nstic  limit.  However,  at  a  water  content  of  7?i  to  8%,  decreasing 
density  had  an  effect  of  increasing  erosion  by  over  1.5  times.  This  wa¬ 
ter  content  was  below  both  the  modified  optimum  of  10?o  and  the  plastic 
limit  of  15%. 

89.  Data  from  Material  B  was  obtained  from  tests  performed  on  mater¬ 
ial  compacted  at  its  optimum  water  content  of  15%,  which  was  6%  below 
the  plastic  limit  of  21%.  The  density  varied  from  90?i  to  98%  of  stan¬ 
dard  Proctor  density  (97  PCF  to  106  PCF).  For  Material  B  the  percent 
erosion  per  unit  time  decreased  from  0.33  %/min  to  0.08  ?0/min  for  this 
increase  in  density. 

90.  In  evaluating  this  data  it  should  be  noted  that  specimens  used 
for  Lewis  and  Schmidt's  study  were  compacted  by  both  the  standard  and 
modified  Proctor  methods.  Therefore,  the  compaction  energy  may  also  in¬ 
fluence.  erosion  for  fabric  sensitive  materials  such  as  clays. 

91.  Effect  of  Eroding  Fluid  Concentration :  Figure  20  shows  that  in¬ 
creasing  the  ionic  concentration  of  the  eroding  fluid  had  little  effect 
on  the  erodibility  of  the  clay  core  material  (Material  B),  but  decreased 
the  erodibility  of  the  silt  core  materials  (Materials  D  and  E). 

92.  Field  observations  and  other  laboratory  investigations  have  also 
indicated  that  a  marked  decrease  in  the  amount  of  erosion  is  expected 
for  an  increase  in  the  ionic  concentration  of  the  eroding  fluid.  Arula- 
nandan,  et.  al.  (1975)  provided  an  explanation  for  this  behavior.  They 
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noted  that  when  the  ionic  concentration  of  the  fluid  is  lower  than  the 


ionic  concentration  of  the  pore  fluid,  water  may  move  into  the  surface 
of  the  soil  by  osmosis.  This  increase  in  water  content  would  increase 
the  interparticle  distance  at  the  surface,  and  thus  reduce  bonding  be¬ 
tween  particles  such  that  the  erodibility  of  the  material  would  be  in¬ 
creased.  This  explanation  is  in  accord  with  the  explanation  offered 
previously  for  the  increase  in  erosion  observed  when  the  molding  water 
content  is  above  the  optimum  water  content  of  the  material. 

93.  On  the  other  hand,  similar  investigations  by  Arulanandan,  et. 
al.  (1973)  have  shown  that  although  two  eroding  fluids  may  have  the  same 
ionic  concentration,  calcium  cations  have  a  larger  effect  on  reducing 
the.  amount  of  erosion  than  sodium  cations.  Moreover,  it  is  interesting 
to  note  that  although  the  total  concentrations  of  the  two  eroding  fluids 
used  in  this  study  were  different  (3  and  6  meq/1),  the  amount  of  calcium 
cations  present  in  each  eroding  fluid  was  the  same  (2  meq/1).  There¬ 
fore,  it  is  possible  that  the  relative  erosion  insensitivity  of  the  ma¬ 
terials  tested  to  different  ionic  concentrations  was  due  to  the  eouiva- 
lent  amount  of  calcium  cations  present  in  the  eroding  fluids. 

94.  Effect  of  Using  Standard  Concrete  Sands  as  Filter  Material :  In 
1977,  Sherman  proposed  the  use  of  standard  concrete  sand  (ASTM  C- 33 )  as 
a  filter  for  all  fine-grained  materials  (Wilson  and  Marsal,  1979).  Re¬ 
sults  of  the  limited  amount  of  data  from  the  filter  tests  performed  for 
this  report  suggested  that  standard  concrete  sand  filters  are  effective 
in  preventing  erosion  of  a  highly  erosive  silt  (Material  E)  at  pressure 
heads  of  up  to  40  meters.  Therefore,  Sherman's  proposition  may  repre- 
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sent  a  simple  and  more  economical  alternative  for  the  protection  of 
cracked  embankment  dam  core  materials  than  present  filter  design  cri¬ 
teria  based  solely  upon  the  grain  size  distribution  of  intact  core  and 
transition  materials  (Table  3). 

95.  Effect  of  Assumed  Erosion  Channel  Shape:  Appendix  D  presents 
the  results  of  a  numerical  analysis  to  investigate  the  validity  of  as¬ 
suming  that  the  triaxial  erosion  test  flow  channel  eroded  as  a  constant 
width  rectangle.  It  may  be  remembered  that  it  was  assumed  that  only  the 
thickness  of  the  rectangular  flow  channel  varied  with  time  to  account 
for  the  soil  volume  being  eroded.  The  results  of  calculations  performed 
using  this  simple  model  were  numerically  compared  to  results  obtained 
using  two  different  flow  channel  shape  models:  1)  a  constant  width  el¬ 
liptical  shape,  and  2j  a  variable  width  elliptical  shape  having  both  the 
width  and  thickness  varying  with  time  to  account  for  the  soil  volume  be¬ 
ing  eroded. 

96.  The  results  of  the  numerical  comparison  showed  that  parameters 
calculated  using  the  simple  constant  width  rectangle  shape  model  varied 
by  less  than  ±1.>*„  from  the  parameters  calculated  using  either  of  the  two 
more  sophisticated  shape  models.  This  variation  was  felt  to  be  insig¬ 
nificant  for  the  results  presented  in  this  report.  Therefore,  the  sim¬ 
ple  rectangular  channel  was  used. 

97.  Effect  of  Assum  *  Hydraulic  Behavior :  Appendix  E  presents  the 
results  of  an  experimental  study  to  observe  and  characterize  the  flow 
occuring  through  a  plexiglas  model  of  the  triaxial  erosion  test  flow 
channel.  Observations  of  these  hvd:aulic  experiments  indicated  that 
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flow  separation  could  be  introduced  during  an  erosion  test  at  the  top  of 
the  erosional  channel  at  the  abrupt  transition  between  the  narrow  in¬ 
jecting  flow  nozzle  and  the  rapidly  enlarging  erosional  channel.  This 
flow  separation  was  not  severe  at  gradients  below  approximately  35  and 
for  crack  thicknesses  less  than  about  0.7  cm. 

98.  An  equation  for  duct  entrance  effects  on  fluid  wall  shear  stress 
was  used  to  estimate  a  30°o  increase  over  laminar  wall  shear  due  to  flow 
separation  effects  (Olson.  1980).  This  30%  increase  in  shear  was  ac¬ 
counted  for  in  the  flow  parameter  ^  in  Eq.5  of  PART  V. 

99.  Some  additional  experiments  were  performed  to  investigate  opti¬ 
mization  of  the  triaxial  erosion  test  hydraulic  design  so  that  flow  sep¬ 
aration  would  be  lessened.  Replacing  the  top  platen  and  injecting  flow 
nozzle  with  a  top  platen  having  a  gradual  transition  conical  flow  en¬ 
trance  resulted  in  a  marked  suppression  of  flow  separation.  Some  addi¬ 
tional  study  is  recommended  to  optimize  the  hydraulic  design  of  the  tri¬ 
axial  erosion  test  apparatus. 
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PART  IX;  SUMMARY  AND  CONCLUSIONS 


A  triaxial  erosion  test  was  used  to  study  the  influence  of  the  engi¬ 
neering  properties  of  soils  and  chemistry  of  the  eroding  fluid  on  the 
erodibility  of  typical  embankment  dam  core  materials.  The  test  modeled 
the  flow  of  water  through  a  cracked  element  of  soil  under  a  state  of 
stress  and  hydraulic  conditions  similar  to  those  existing  in  the  core  of 
a  dam.  The  following  conclusions  were  reached  based  on  the  results  of 
the  investigation: 

1)  High  erosion  took  place  in  silt  soils,  moderate  erosion  took  place 
in  clayey  silt  soils  and  effectively  no  erosion  took  place  in  clay 
soils.  These  results  are  in  agreement  with  results  obtained  by  others 
using  the  pinhole,  test,  and  other  similar  test  procedures. 

2)  The  erodibility  of  silt  materials  is  greatly  influenced  by  varia¬ 
tion  in  the  molding  water  content.  However,  the  effect  of  molding  water 
content  on  erodibility  is  less  pronounced  for  clayey  materials.  These 
results  suggest  that  control  of  the  compaction  water  content  during  con¬ 
struction  may  reduce  the  erosion  potential  of  embankments  with  silt 
cores,  but  not  for  embankments  with  ciay  cores. 

3)  The  water  content  at  the  plastic  limit  has  been  suggested  in  the 
literature  as  an  indicator  of  the  moisture  content  at  which  minimum  ero¬ 
sion  will  occur.  In  contrast,  the  results  of  this  investigation  indi¬ 
cate  that  the  optimum  compaction  water  content  is  a  better  indicator 
than  the  plastic  limit  of  the  moisture  content  at  which  minimal  erosion 
will  occur.  This  inconsistancv  can  be  explained  when  it  is  noted  that 
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for  some  of  the  materials  cited  In  the  literature,  the  water  content  at 
the  plastic  limit  was  very  close  to  the  optimum  compaction  moisture  con¬ 
tent  . 

4)  Limited  data  obtained  indicated  that  density  has  little  effect  on 
the  erosion  potential  of  silt  embankment  dam  core  material  compacted  at 
the  optimum  molding  water  content.  However,  clayey  materials  were  some¬ 
what  sensitive  to  variations  in  density;  for  clays  erodibility  increased 
as  density  decreased. 

5)  Tests  also  showed  that  the  ionic  concentration  measurably  influ¬ 
enced  the  erodibility  of  silt  core  materials  but  had  less  effect  on  the 
erodibility  of  clay  core  materials.  The  higher  the  ionic  concentration 
of  the  eroding  fluid,  the  lower  the  erosion. 

6)  Tests  indicated  that  standard  concrete  sand  material  (ASTM  C-33) 
as  a  downstream  filter,  clogs  core  cracks  by  collecting  core  material 
particles  downstream.  Such  downstream  filters  effectively  prevented  ero¬ 
sion  even  for  highly  erosive  materials  under  hydraulic  heads  of  up  to  40 
meters.  Such  sand  filter  materials  may  provide  more  economical  protec¬ 
tion  of  cracked  embankment  dam  core  materials  than  filters  designed  us¬ 
ing  criteria  based  upon  the  grain-size  distribution  of  core  and  tran¬ 
sition  materials. 

7)  It  was  found  that  the  erosion  rate,  i,  asymptotically  approached 
a  maximum  value  for  increasing  values  of  fluid  shear  stress.  If  a  maxi¬ 
mum  value  of  erosion  rate  does  in  fact  exist  as  a  material  constant  un¬ 
der  known  compaction  and  enviromental  conditions,  it  would  be  a  useful 
material  parameter  for  embankment  dam  design. 


8)  It  is  recommended  that  more  research  be  conducted  to: 

a)  Optimize  the  hydraulic  characteristics  of  the  triaxial  erosion 
test  apparatus. 

b)  Further  study  the  inter-relationship  between  moisture  content, 
density  and  the  erosion  rate. 

c)  Further  study  the  possibility  that  a  maximum  erosion  rate  ex¬ 
ists,  and  how  this  concept  may  be  used  in  embankment  dam  design. 
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erosion  test  program  parameters 
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SUMMARY  SHOWING  THE  ERODIBILITY  OF  MATERIALS 
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Fig.  2  Photographs  of  Erosion  Test  Equipment  and  Erosion 
Test  Cell. 


. 


Fig.  3  Schematic  Representation  of  Triaxial  Erosion  Test  Apparatus 


Fig.  5  Photographs  of  a)  mold  used  to  prepare  reconstituted 
specimens  showing  the  nozzle  insert  and  b)  specimen 
extruder  modified  to  form  a  crack  in  the  erosion 
specimen. 
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^ig.  6  Standard  Proctor  compaction  curve  and  gram 
size  distribution  for  Material  A. 
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size  distribution  for  Material  D. 
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Fig.  10  Standard  Proctor  compaction  curve  and  grain 
size  distribution  for  Material  E. 

(Teton  Dam  Embankment  Material) 
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APPENDIX  C:  PROCEDURE  FOR  OBTAINING  SOIL  WATER  EXTRACTS 


This  appendix  describes  the  procedure  recommended  by  Sherard  (1972) 
and  extensively  described  by  Richards  (1954)  for  obtaining  soil-water 
extracts.  These  extracts  are  used  to  determine  the  amount  of  metallic 
cations  (sodium,  calcium,  magnesium  and  potassium)  present  in  the  soil, 
expressed  in  milliequivalents  per  liter.  These  ionic  concentrations  may 
be  used  to  determine  the  relative  activity  of  sodium  ions  in  exchange 
reactions  with  soil. 

A  high  activity  of  sodium  cations  relative  to  the  other  ions  in  the 
porewater  has  been  found  to  correlate  to  the  incidence  of  erosion  damage 
in  many  embankment  dams.  One  widely  used  relationship  to  quantify  this 
activity  is  the  Sodium  Adsorption  Ratio  (SAR),  defined  as: 


Na  (1) 

SAR  =  -  -  — 

/  0.5  (  Ca  +  Mg  ) 

The  SAR  may  be  plotted  versus  the  total  cation  concentration  of  the 
eroding  fluid  to  possibly  predict  the  erosion  failure  of  a  dam  embank¬ 
ment.  To  obtain  the  soil-water  extracts  required  to  compute  the  SAR, 
the  following  procedure  is  recommended: 

1.  Distilled  water  is  added  to  a  250  gram  soil  sample  and  mixed  with  a 
spatula  until  the  consistency  of  the  soil  paste  is  about  the  same  as  at 
the  liquid  limit.  The  mixing  process  may  be  facilitated  if  the  soil  is 
first  allowed  to  air  dry  and  then  passed  through  a  U.S.  Standard  No.  70 
(2mm)  sieve. 
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2.  The  soil  paste  is  allowed  to  sit  for  a  few  hours  in  order  to  allow 
equilibrium  to  occur  between  the  salts  in  the  pore  water  and  the  salts 
on  the  cation  complex. 

3.  After  allowing  the  paste  to  sit  for  a  few  hours,  the  paste  should  be 
checked  to  see  that  free  water  does  not  collect  on  the  soil  surface.  If 
the  paste  is  too  wet,  dry  soil  may  be  added.  The  paste  should  not 
stiffen  or  loose  its  glister.  If  the  paste  is  too  dry,  it  may  be  remix¬ 
ed  with  additional  water.  The  consistency  of  the  paste  should  be  about 
the  same  as  at  the  liquid  limit. 

A.  The  soil  paste  is  transfered  to  a  vacuum  filter  funnel  system,  as 
shown  in  Figure  Cl.  Vacuum  is  applied  and  a  small  quantity  (10  to  25 
milliliters)  of  soil-water  extract  is  collected  in  a  bottle  or  test 
tube.  The  collection  process  is  stopped  when  air  begins  to  pass  through 
the  filter. 

5.  The  amount  of  metallic  cations  is  determined  by  chemical  analysis. 
A  detailed  description  of  the  procedures  available  for  determining  pore 
water  chemistry  and  other  properties  of  the  soils  may  be  found  in  texts 
edited  by  Richards  (1954)  and  Black  (1965). 
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APPENDIX  D:  NUMERICAL  ANALYSIS  OF  THE  EFFECT  OF  ASSUMED  FLOW 


CHANNEL  SHAPE  ON  TRIAXIAL  EROSION  TEST  RESULTS 

Studies  were  performed  to  evaluate  the  effect  of  flow  channel  shape 
on  parameters  calculated  for  the  triaxial  erosion  test.  As  discussed  in 
the  report,  test  results  were  computed  assuming  the  flow  channel  shape 
was  rectangular  during  a  test..  This  rectangular  channel  was  assumed  to 
have  a  constant  width;  only  the  thickness  varied  with  time  to  account 
for  the  soil  volume  being  eroded.  Possibly  more  accurate  assumptions 
could  also  have  been  made  by  assuming  the  flow  channel  eroded  as  an  el* 
lipse  with  either  a  constant  width  or  with  both  width  and  thickness  var¬ 
ying  with  time.  Therefore  comparison  calculations  were  performed  by  as¬ 
suming  these  differently  shaped  flow  channels.  Flow  channel  size 
parameters  compared  were:  width  (w),  thickness  (d),  hydraulic  diameter 

(D.  ),  and  erosional  surface  area  (A  (t)).  The  effects  of  the  flow 
n  s 

channel  shape  assumption  on  the  erosion  rate  and  fluid  shear  stress  were 
also  compared. 

Method  of  Study 

Erosion  test  values  calculated  for  the  two  different  elliptical  shape 
flow  channels  were  compared  to  test  values  calculated  for  the  constant 
width  rectangular  shape  flow -channel  on  the  basis  of  equal  cross-sec¬ 
tional  areas,  (t).  This  was  done  because  the  area  of  any  cross- 

section  is  independent  of  cross-sectional  shape. 

The  constant  width  solutions  were  selected  based  on  observation  of 
the  eroded  flow  channels  after  a  test.  These  observations  indicated 
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that  the  increase  in  crack  width  was  much  smaller  than  the  increase  in 


crack  thickness.  The  variable  width  solution  attempted  to  model  this  by 
using  an  empirical  function  to  generate  values  of  the  small  increase  in 
crack  width  during  an  erosion  test. 


Constant  Width  Rectangular  Shape  Flow  Channel :  The  solution  for  the 

constant  width  rectangular  shape  flow  channel  is  described  in  PART  IV  of 
this  report.  The  expressions  for  the  channel  thickness  d  and  the  chan¬ 
nel  area  A  (t)  as  functions  of  time  were: 
s 


d  (t)  =  Ax  (t)/w  (Dl) 

and 

A  (t)  ■  1  *  2  *  (  A  (t)  /w  +  w  )  (D2) 

s  e  x 

where: 

lp  *  eroded  length  of  specimen,  in  cm, 

d  (t)  *  crack  thickness  at  time  *  t,  in  cm, 
w  =  original  crack  width  at  time  *  0,  in  cm. 

A&  (t)  *  crack  surface  area  at  time  *  t,  in  sq.  cm 

A^  (t)  =  crack  cross-sectional  area  at  time  *  t,  in  sq.  cm. 
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Constant  Width  Elliptical  Shape  Flow  Channel:  The  constant  width  el¬ 


lipse  solution  was  based  on  mensuration  formulae  which  were  used  to  ob¬ 
tain  values  of  d  and  A  (t)  as  functions  of  w  and  A  (t).  These 

s  x 

mensuration  formulae  were: 


d  (t)  =  (  4/it  )  *  Ax  (t)/w 


(D3) 


and 


A  (t)  =  1  *  (  it/2  )  *  (  w+d  )  *  (  64-  3*  Z*  )  /  (  64-  16*  Zl  )  (D4) 

s  e 


where: 


Z  «  [  w-  d  (t)  ]  /  l  w+  d(t)  ] 


Variable  Width  Elliptical  Shape  Flow  Channel :  This  solution  utilized 
basically  the  same  mensuration  formulae  as  above,  but  required  an  addi¬ 
tional  function  to  generate  values  of  w  and  d  as  the  cross-sectional 
area  A^  (t)  varied  with  time.  For  our  computations,  we  empirically 

derived  a  function  based  on  observed  behavior  of  erosion  channels  such 
that: 


Aw  »  l  AAx  (t)/At  j  *  Ad  *  RATE 


(D5) 


where: 


Aw  *  incremental  change  in  crack  width  in  cm. 

Ad  *  incremental  change  in  crack  thickness  in  cm, 

I* 


AAx  /At  =  change  in  cross-sectional  area  with  time,  cm2/min, 
and 

RATE  =  constant,  min/cm*. 

[  The  value  of  AAx  (t)/At  was  equal  to  P  (the  slope  of  the 

cumulative  weight  eroded  vs.  time  curve)  divided  by  p  (the 
specimen  dry  density)  and  lg  (the  eroded  length  of  the  spec¬ 
imen  ) .  ] 

Equation  D5,  with  a  suitable  choice  for  the  value  of  RATE,  could  fit 
a  wide  range  of  crack  dimensions  measured  at  the  end  of  an  erosion  test. 
This  equation  was  incrementally  solved  by  using  consecutive  values  of 
Ax  (t)  and  by  starting  with  initial  values  of  w  and  d  at  time  =  0. 

Once  these  incremental  values  were  known,  the  elliptical  mensuration 
formulae  could  be  used  to  compute  flow  channel  size  parameters. 

The  RATE  parameter  was  used  to  evaluate  the  time  history  of  the  the 
shape  of  the  flow  channel  based  on  measurements  of  flow  channel  dimen¬ 
sions  at  the  end  of  a  triaxial  erosion  test.  Therefore  values  for  RATE 
necessary  to  match  the  the  test  results  were  calculated  analytically  us¬ 
ing  measured  test  end  values  of  d,  w,  time  (or  Ax  (t))  and  AA^  /At. 

An  example  of  this  tabulation  is  given  in  Table  Dl.  Tables  D2  through 
D4  present  the  measured  test  end  crack  dimensions  and  best  fit  values  of 
RATE  for  most  of  the  specimens  tested. 
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Results  of  Study 


Constant  Width  Ellipse:  A  numerical  comparison  of  erosion  test  val¬ 
ues  calculated  using  the  constant  width  ellipse  solution  and  the  con¬ 
stant  width  rectangle  solution  is  presented  in  Tables  D5  through  D7  for 
each  of  the  materials  studied.  The  analysis  indicated  that  for  all  of 
the  materials,  the  constant  width  ellipse  crack  size  parameters  vary 
from  the  rectangular  by  less  than  ±  15%.  It  should  also  be  noted  that 
the  ratio  of  the  results  between  solutions  were  about  the  same  for  all 
of  the  materials  studied. 

Figure  D1  shows  graphically  the  effect  of  the  assumed  flow  channel 
shape  on  values  of  erosion  rate,  t,  versus  fluid  shear  stress,  t.  Ero¬ 
sion  rate  values  for  the  constant  width  ellipse  solution  were  consis¬ 
tently  higher  than  the  values  for  rectangular  solution  of  the  materials 
tested.  Plots  of  erosion  rate  verses  fluid  shear  stress  for  Material  B 
were  not  sensitive  to  the  choice  of  assumed  flow  channel  shape  because 
the  measured  erosion  rates  were  small. 

Erosion  test  values  for  the  highly  erosive  Material  E  were  more  sen¬ 
sitive  to  the  choice  of  assumed  flow  channel  shape  for  the  less  erosive 
materials.  The  slopes  of  the  constant  width  ellipse  solution  were  sig¬ 
nificantly  greater  than  the  slopes  of  the  rectangular  solution  for  this 
material.  Calculated  erosion  values  for  the  two  solutions  diverged  mod¬ 
erately  with  increasing  values  of  fluid  shear  stresses. 

Calculated  erosion  values  for  the  two  solutions  fo*-  the  moderately 
erosive  Material  D  were  generally  parallel  except  at  low  fluid  shear 
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stresses.  At  low  shear  stress  the  slope  of  the  constant  width  ellipse 
solution  was  slightly  greater  than  the  slope  of  the  rectangular  solu¬ 
tion. 

Variable  Width  Ellipse :  The  numerical  comparison  of  the  variable 
width  ellipse  solution  to  the  constant  width  rectangle  solution  is  pre¬ 
sented  in  Tables  D8  through  DIO  for  each  of  the  materials  studied.  The 
variable  width  ellipse  solution  also  varied  from  the  rectangle  solution 
by  less  than  ±  15%  for  all  of  the  materials  tested. 

Figure  D1  presents  graphical  comparisons  of  the  effect  of  assumed 
flow  channel  shape  on  values  of  erosion  rate,  i,  versus  fluid  shear 
stress,  t .  For  each  of  the  tested  materials  the  variable  width  ellipse 
solution  resulted  in  curves  higher  than,  but  generally  parallel  to,  the 
curves  for  the  rectangle  solution.  Further,  it  may  be  seen  in  Figure  D1 
that  the  shapes  of  the  curves  for  the  two  solutions  are  quite  similar. 
Conclusions 

The  results  of  the  numerical  analyses  indicated  that  the  simple  con¬ 
stant  width  rectangular  solution  was  reasonable.  Another  encouraging 
finding  is  that  erosion  values  calculated  using  the  possibly  most  accu¬ 
rate  3hape  assumption  of  a  variable  width  ellipse  may  be  related  to  the 
constant  width  rectangle  by  using  a  simple  constant. 

An  alternative  test  method  that  would  avoid  all  of  these  shape  as¬ 
sumptions  would  be  to  start  the  erosion  test  with  a  circular  flow  chan¬ 
nel  shape.  Such  tests  were  actually  performed  by  making  flow  channels 
in  specimens  by  drilling  with  a  small  (0.7  cm  diameter)  auger  bit.  Ob- 
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servations  indicated  that  the  resulting  flow  channel  eroded  circularly. 
The  mathematics  for  calculating  the  change  in  size  of  an  eroding  circu¬ 
lar  flow  channel  are  very  simple  and  no  flow  channel  shape  assumptions 
are  required.  Also,  since  the  aspect  ratio  of  the  circular  channel  does 
not  vary,  computation  of  fluid  wall  shear  stress  as  a  function  of  aspect 
ratio  is  simplified.  Additional  study  is  required  to  see  if  erosion  pa¬ 
rameters  measured  using  a  circular  flow  channel  compare  favorably  with 
values  measured  using  a  rectangular  flow  channel. 


Table  01.  Example  Tabulation  for  Determining  the  Value  of  RATE 
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Table  D-2.  Tabulation  of  Test  End  Dimensions  for  Material  B  Specimens. 
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Table  D-5.  Tabulation  of  Numerical  Comparison  Between  Constant  Width 

Ellipse  and  Constant  Width  Rectangle  Solutions  for  Specimen  B-117. 
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Table  D-6.  Tabulation  of  Numerical  Comparison  Between  Constant  Width  Ellipse 
and  Constant  Width  Rectangle  Solutions  for  Specimen  D-175 . 


i 

i 

i 

1 

1 

H*  1 

i 

m  i 

• 

<  i 

PO 

• 

S  1 

r— 

i 

U  | 

• 

• 

• 

i 

(A  l 

f— 

Z  i 

© 

<  l 

o  1 

►-  m 

l 

1 

X 

1 

O  1 

X  Ui 

1 

O  1 

< 

l 

jy 

X  l 

<0 

® 

*—  1 

sz 

C  l 

<0 

<0 

« 

*»  1 

o  < 

\  i 

• 

• 

• 

1 

XX 

X  1 

© 

© 

© 

<  < 

C  i 

1 

*-  CL 

1 

“  1 

u 

1 

^  1 

Ui  -i 

1 

o  » 

X  < 

1 

Z  1 

o 

1 

© 

s 

Ui  1 

u.  — 

:*  1 

c 

© 

o 

_J  1 

QH. 

N  1 

• 

• 

• 

1 

CL 

*  • 

C  l 

©  - 

l 

Ui  • 

-  -J 

1 

Q  • 

N»  -i 

1 

O  l 

<Ui 

1 

X  • 

x 

1 

U  1 

1 

Os 

Cs 

Os 

1 

e  i 

r*» 

>• 

1 

S  1 

• 

• 

• 

1 

C  l 

© 

© 

© 

1 

o  < 

1 

u  • 

1 

1 n  % 

N  • 

1 

uj  1 

©  « 

1 

a.  i 

1 

1 

<  i 

*»  1 

1 

PP*  1 

Ui 

1 

<7?  i 

1 

c  u 

Z  1 

i 

Ui  <  < 

-  o  i 

CM 

CM 

a  • 

1 

©  u.  Li  <A  *l 

• 

• 

• 

U  1 

ui 

1 

o  z  z  <  a  i 

*— 

P*. 

X  • 

(A 

Z  9< 

(A  l 

>0 

CO 

<A©  1 

a. 

• 

w  m 

ui  (A 

< 

—  m  • 

>  1 

XX 

1 

(/)</}  • 

-J 

H  ( 

<  Ui 

1 

>  <  • 

_  • 

“  ^ 

C 

I 

-i  i 

©Ui 

<A  1 

(A  Ui 

z  < 

ZZ  1 

4A 

CM 

s 

<  U.  1 

Z  1 

z 

o- 

©O  i 

• 

• 

ZQ  I 

w  x 

Ui  1 

-i  < 

>© 

1 

*“ 

CM 

CM 

<  • 

O  l 

<  X 

x 

1 

in  % 

zo 

• 

u  t 

1 

u.  z  i 

w  — 

>  i 

•  CL 

1 

Ow  i 

Z  5 

X  1 

• 

1 

C  1 

a.  x 

«• 

Z  1 

PO 

tA  -j  i 

3  W 

1 

-q 

9 

3CO  1 

« 

• 

• 

*■»  z  • 

UiZ 

J  • 

«• 

1 

CM 

CM 

-J<  1 

fl-< 

■>  1 

JZ 

1 

3K  1 

</>  — 

O  1 

Ui  U 

1 

<A<  • 

(A 

<A  1 

1 

im  a  < 

z 

1 

* 

1 

Z  I 

o 

l 

u 

^  Z  ■ 

CM 

W  1 

o 

• 

— 

CO  i 

• 

• 

• 

X  1 

l 

x 

I 

• 

CM 

N 

<  • 

Z  l 

i 

CL  1 

—  1 

i 

Z  1 

Z  1 

i 

©  i 

S  1 

Ui 

i 

u  • 

O  1 

cu 

Z  1 

1 

1 

Ui  <  < 

-O  * 

>0 

Q  1 

c  • 

a  uj  <a  *  i 

• 

• 

• 

• 

Ui 

©zz 

Z  9< 

“sr ; 

92 

O 

Cs 

>0  l 

Z  (A 

—  (A 

i 

<M  l 

<  Z 

i 

l 

5- 

o 

i 

i 

1 

Z  l 

Ui 

cs  z 

l< 

llll 

CM 

OS 

* 

©  • 

<  < 

>© 

1 

fJ 

JZ 

Z 

1 

(A  l 

9  < 

1 

O  i 

OO. 

1 

X  1 

z 

x 

1 

ui  • 

1 

<  * 
u  q 

3 

W 

CM 

us 

U  1 

o  < 

1 

CM* 

CM 

(Si 

:  i 

Ui  X 

5 

< 

3  • 

X  u 

1 

w  • 

* 

1 

9  • 

u 

X  1 

Os 

• 

«• 

CU  i 

• 

• 

• 

u»  i 

X 

i 

C 

• 

C  i 

•• 

i 

i 

W  1 

i 

(A 

i 

<  » 

<A< 

i 

s  I 

i 

1 

CZ 

z  • 

1 

u< 

-  U  1 

CNi 

Cs 

*Ts 

1 

X  •  • 

• 

* 

• 

1 

*  • 

<  ©  1 

CM 

1 

y  *■ 

;A  i 

1 

<  u 

I 

1 

Z  - 

I 

1 

O'-A 

I 

D15 


Table  D-7.  Tabulation  of  Numerical  Comparison  Between  Constant  Width  Ellipse 
and  Constant  Width  Rectangle  Solutions  For  Specimen  E-102. 
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Table  D-8.  Tabulation  of  Numerical  Comparison  Between  Variable  Width  Ellipse 
and  Constant  Width  Rectangle  Solutions  For  Specimen  B-117. 
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APPENDIX  E:  EXPERIMENTAL  STUDY  OF  TRIAXIAL  EROSION  TEST 


HYDRAULIC  BEHAVIOR 

A  series  of  hydraulic  experiments  was  performed  to  observe  flow 
through  a  model  of  the  erosional  flow  channel.  These  experiments  were 
performed  to  investigate  the  validity  of  the  assumption  of  laminar  flow 
occuring  in  the  triaxial  erosion  test.  Additionally  some  experiments 
were  performed  to  possibly  optimize  the  hydraulic  design  of  the  top  and 
bottom  platens  presently  used  in  the  triaxial  erosion  test. 

Experimental  Apparatus 

A  clear  rectangular  plexiglas  channel  was  constructed  to  model  the 
erosional  flow  channel.  The  model  channel  dimensions  were  11.5  cm  long, 
by  2.9  cm  wide,  by  1.2  cm  thick.  These  dimensions  approximately  corre¬ 
spond  to  the  values  of  1  ,  w  and  d  used  for  the  tests  performed  in 

this  report.  The  model  channel  thickness  could  be  varied  by  placing 
plexiglas  inserts  0.5  cm  thick  into  the  channel.  Gradients  of  approxi¬ 
mately  10,  35,  50  and  100  m/m  were  applied  to  each  of  the  three  model 
channel  configurations  possible:  0.2  cm,  0.7  cm  and  1.2  cm  thick  by  2.9 
cm  wide. 

The  resulting  flow  regime  was  observed  by  mixing  light  reflective 
polyethylene  beads  with  the  water  flowing  through  the  model  channel.  The 
average  diameter  of  the  beads  was  between  0.03  and  0.05  cm.  The  beads 
had  a  specific  gravity  of  slightly  greater  than  1.  Observations  were 
recorded  by  sketching  and  by  photography. 
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Flows  were  channeled  through  the  plexiglas  model  with  either  the  same 
top  and  bottom  platens  used  for  the  triaxial  erosion  test  or  with  alter¬ 
native  platens.  The  platens  were  sealed  to  the  model  channel  using  mod¬ 
eling  clay.  Figure  El  shows  a  photograph  of  the  experimental  apparatus 
during  use. 

To  see  the  effect  of  entrance  and  exit  geometry  on  the  flow  regime,  a 
platen  with  a  gradual  transition  conical  flow  entrance  (0.8  cm  expanding 
out  to  2.2  cm  in  diameter  over  a  length  of  1.2  cm)  was  substituted  for 
either  of  the  top  or  bottom  platens.  Figure  E2  shows  details  of  the 
conical  flow  entrance  plater..  An  "infinitely  large"  exit  was  also  stud¬ 
ied  by  discharging  the  flow  from  the  model  channel  dire  ,.s,1y  into  a  large 
beaker.  This  wa.'  done  for  both  the  original  top  platen  and  the  conical 
flow  entrance  top  platen. 

Experimental  Observations ,  Triaxial  Erosion  Test  Model 

Crack  Thickness  ■  0.2  cm  (This  thickness  is  analogous  to  the  specimen 
flow  channel  thickness  at  the  start  of  the  erosion  test).  Flows  were 
observed  to  be  quite  linear  and  regular  at  all  applied  gradients.  Some 
small  vortices,  about  1  cm  in  diameter,  were  noticed  along  the  edges  of 
the  model  channel  just  below  the  injecting  flow  nozzle. 

Crack  Thickness  f  0.7  cn  (This  approximately  corresponds  to  Ax  (t) 

=2.0  cm*  ).  Two  vortices  were  observed:  one  large  (1.5  cm  wide  by  3 
cm  long)  well-defined  vortex  occured  just  below  one  corner  of  the  in¬ 
jecting  flow  nozzle,  and  a  smaller  (1  cm  wide  by  2  cm  long)  poorly 
formed  vortex  occured  below  and  to  the  other  side  of  the  model  channel 
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opposite  to  the  first  vortex.  This  asymmetric  arrangement  became  more 
and  more  symmetrical  as  applied  gradients  were  increased.  Random  flow 
patterns  with  a  length  of  about  2  cm  occured  just  below  the  vortices. 
The  vortex  fields  took  place  over  about  25™  of  the  area  of  the  model 
channel . 

Crack  Thickness  ®  1.2  cm  (This  approximately  corresponds  to  (t) 

-3.5  cmJ  ).  An  intense  vortex  field  filled  the  top  1/4  of  the  model 
channel.  At  gradients  below  35,  the  flow  in  the  rest  of  the  model  chan¬ 
nel  was  about  equally  divided  between  an  upper  random  non-linear  flow 
region  and  a  lower  linear  flow  region.  As  gradients  were  increased,  the 
vortex  field  in  the  top  of  the  model  channel  became  more  intense,  but 
did  not  change  appreciably  in  extent;  the  random  flow  zone  increased  in 
length  with  a  corresponding  decrease  in  the  linear  flow  zone.  At  a  gra¬ 
dient  of  abcut  100,  the  linear  flow  zone  disappeared  entirely. 

The  observation  of  the  intense  vortex  field  at  the  top  of  the  model 
channel  corresponded  with  test  observations  that  most  erosion  test  spec¬ 
imens  failed  at  a  locally  severely  eroded  area  generally  within  the  top 
1/5  of  the  specimen.  The  intense  vortex  field  may  have  caused  increased 
erosion  due  to  locally  increased  flow  velocities  in  this  region. 

Saturation  of  Flow  Channel  Another  experimental  observation  was  that 
air  pockets  tended  to  form  in  the  plexiglas  model  of  the  flow  channel. 
Such  air  pockets  could  prevent  the  eroding  fluid  from  acting  on  the  ero- 
sional  surface  area  of  a  specimen.  An  effective  solution  to  this  prob¬ 
lem  was  to  pre-fill  the  model  channel  with  fluid  immediately  before  a 
flow  test.  This  was  done  by  1)  holding  the  end  of  the  outflow  line 
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above  the  level  of  the  top  of  the  plexiglas  model,  2)  by  filling  the 
outflow  line  with  fluid,  and  3)  by  allowing  the  trapped  air  to  escape 
from  the  disconnected  inflow  lint.  It  is  recommended  that  the  triaxial 
erosion  testing  procedure  used  in  this  study  be  slightly  modified  to  include 
pre-filling  the  flow  channel  with  the  eroding  fluid  immediately  before  an 
erosion  test. 

Experimental  Observations ,  Optimization  Models 

Gradual  Transition  Conical  Flow  Entrance  Platen  Replacing  the  triaxi¬ 
al  erosion  test  top  platen  and  its  injecting  flow  nozzle  with  a  platen 
with  a  gradual  transition  conical  flow  entrance  had  a  marked  effect  on 
suppressing  the  onset  of  flow  separation  and  minimized  the  formation  of 
vortices.  This  was  probably  because  the  abrupt  transition  between  the 
injecting  flow  nozzle  and  the  model  channel  was  eliminated.  However  it 
should  also  be  noted  that  removing  the  3.5  cm  long  injecting  flow  nozzle 
increased  the  effective  length  of  the  model  channel  from  11.5  cm  to  15 
cm.  This  increase  in  length  may  have  also  tended  to  suppress  the  flow 
separation. 

The  observed  flow  using  the  conical  flow  entrance  was  linear  for 
crack  thicknesses  of  0.2  and  0.7  cm  at  gradients  between  10  and  100. 

Some  small,  1  cm  in  diameter,  vortices  formed  just  below  the  conical 
flow  entrance  at  a  crack  thickness  of  0.7  cm.  Higher  gradients  had  lit¬ 
tle  effect  on  this  flow  regime. 

For  a  crack  thickness  of  1.2  cm,  two  vortices  (3  cm  long  by  1  cm 
wide!  symmetrically  flanked  the  center  line  of  the  model  channel  just 
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beyond  the  conical  flow  entrance.  A  random  flow  area  of  about  3  cm  long 
was  present  below  the  vortices,  followed  by  linear  flow  throughout  the 
rest  of  the  model  channel.  Higher  gradients  had  little  effect  on  this 
regime. 


Alternative  Discharge  Methods:  The  triaxial  erosion  test  bottom  pla¬ 
ten  with  a  0.4  cm  diameter  discharge  was  replaced  either  by  1)  the  coni¬ 
cal  flow  platen  or  by  2)  discharging  directly  into  a  large  beaker.  Lit¬ 
tle  visible  effect  on  the  flow  regime  was  noticed  using  these  different 
discharge  methods.  However,  it  was  noticed  that  coarser  particles  had  a 
tendency  to  lodge  on  the  flat  area  of  the  triaxial  erosion  test  bottom 
platen.  This  did  not  occur  for  the  alternative  bottom  platen  as  this  dead  area 
did  not  exist  (see  Fig.  E2).  To  prevent  the  lodging  of  particles  and  possi¬ 
ble  erosion  channel  sealing,  the  design  of  the  discharge  through  the 
bottom  platen  may  be  modified  in  the  future. 

Numerical  Analysis  of  Hydraulic  Behavior 

The  plexiglas  model,  using  the  triaxial  erosion  test  top  and  bottom 
platens,  revealed  that  flow  separation  increased  as  the  thickness  of  the 
model  channel  increased.  This  flow  separation  would  tend  to  increase 
the  wall  shear  stresses  due  to  locally  increased  flow  velocities.  An 
explicit  derivation  for  the  exact  amount  of  increase  is  not  possible  at 
present.  However,  an  estimate  of  this  increase  in  wall  shear  stress  was 
obtained  by  using  an  equation  for  duct  entrance  effects  on  wall  shear 
(Olson,  1980).  The  results  of  this  analysis  indicated  that  laminar  flu¬ 
id  wall  shear  could  be  increased  up  to  about  30%  due  to  flow  separation. 
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In  our  analysis  of  the  assumed  fully  developed  laminar  flow  through 
the  triaxial  erosion  test  flow  channel,  the  wall  shear  stress,  t,  varied 
with  the  aspect  ratio  of  the  flow  channel  (Olson,  1980).  Thus  the  wall 
shear  stress,  t,  was  equal  to  10.6  Vp/D^  at  the  start  of  the  erosion 

test  when  the  aspect  ratio  was  1/10,  and  decreased  to  a  value  of  8 
Vp/D^  when  the  flow  channel  eroded  to  an  approximately  circular 

shape.  The  average  value  of  wall  shear  stress  t  was  thus  equal  to  9.3 

Vp/D,  .  Therefore  our  test  values  for  T  were  increased  by  30*  to  ac- 
n 

count  for  the  flow  separation  observed  in  the  plexiglas  model  of  the 
triaxial  erosion  test  flow  channel.  This  correction  results  in  the  wall 
shear  stress  t  being  equal  to  12  Vp/D^  .  This  last  value  was  used  in 

the  computation  of  wall  shear  stresses  reported  in  the  text. 

Conclusions 

Observations  indicated  that  as  the  flow  channel  erodes,  flow  separa¬ 
tion  occurs  at  the  abrupt  transition  between  the  top  platen  injecting 
flow  nozzle  and  the  flow  channel.  The  experiments  indicated  that  upper 
bound  solutions  for  wall  shear  stresses  due  to  the  assumed  fully  devel¬ 
oped  laminar  flow  are  more  valid  than  solutions  based  on  average  condi¬ 
tions.  To  reflect  this,  wall  shear  stresses  were  increased  by  30*  over 
average  values  to  obtain  the  results  presented  in  this  report. 

Some  experiments  were  also  conducted  to  investigate  the  triaxial  ero¬ 
sion  test  apparatus  hydraulics  so  that  flow  separation  effects  would  be 
lessened.  These  experiments  indicated  that  by  changing  the  top  platen 
having  an  abrupt  flow  nozzle  to  a  top  platen  with  a  gradual  transition 

conical  flow  entrance,  a  narked  suppression  of  flow  separation  resulted. 


E7 


However,  increasing  the  size  or  changing  the  geometry  of  the  flow 
discharge  through  the  bottom  platen  had  little  visible  effect  on  flow 
separation. 

To  improve  the  accuracy  of  the  hydraulic  analysis  for  the  fluid  wall 
shear  stresses,  further  study  is  needed  to  optimize  the  configuration  of 
the  entrance  to  the  flow  channel.  The  goal  of  this  optimization  would 
be  to  lessen  the  flow  separation  occur ing  in  the  present  entrance  con¬ 
figuration.  Another  goal  would  be  to  prevent  flow  from  occuring  between 
the  specimen  and  its  membrane.  The  present  top  platen  and  injecting 
flow  nozzle  do  prevent  flow  from  occuring  between  the  specimen  and  its 
membrane,  but  could  better  control  flow  separation  in  the  transition  be¬ 
tween  the  flow  entrance  and  the  flow  channel. 

It  is  recommended,  as  previously  discussed,  that  the  triaxial  erosion 

j 

testing  procedure  be  slightly  modified  to  include  pre-filling  the  flow 
channel  with  the  eroding  fluid  immediately  before  an  erosion  test. 


Fig.  El  Experimental  Apparatus  Used  to  Investigate  Triaxial 
Erosion  Test  Hydraulic  Behavior • 
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Kin 


appendix  f 

estimates  of  field  erosion  rates 
for  THE  TETON  DAM  FAILURE  AND 
POSSIBLE  DESIGN  APPLICATIONS 


FI 


APPENDIX  F:  ESTIMATES 
OF  FIELD  EROSION  RATES 
FOR  THE  TETON  DAM  FAILURE  AND  POSSIBLE 
DESIGN  APPLICATIONS 


Background 

Based  on  field  observations  during  the  failure  of  the  Teton  Dam, 
calculations  presented  in  this  Appendix  were  made  to  obtain  estimates 
of  the  erosion  rate  which  may  have  occured  while  an  erosion  channel  was 
enlarging  during  the  dam  failure.  These  estimates  were  compared  to 
values  of  the  maximum  erosion  rate  obtained  using  the  Triaxial  Erosion 
Test  to  observe  if  laboratory  data  corresponded  with  field  observations. 

The  following  paragraphs  summarize  key  erosion  events  presented  in  a 
report  on  the  failure  of  the  Teton  Dam  on  the  morning  of  June  5,  1976 
(U.S.  Department  of  the  Interior,  1977). 

Earliest  Time  at  which  Significant  Erosion  may  have  begun: 

Testimony,  presented  on  page  F-37  of  the  aforementioned  report,  established 
that  no  serious  leaks  or  seepage  were  observed  during  a  routine  inspection 
at  9:00  p.m.  on  June  U ,  1976.  Therefore,  the  earliest  time  at  which 
significant  erosion  could  have  begun  would  be  subsequent  to  9:00  p.m. 
on  June  U,  1976  (Page  35). 

Latest  Time  at  which  Significant  Erosion  may  have  begun:  The 
first  leak  was  observed  at  about  7:^5  a.m.,  on  June  5,  1976.  This  leak 
occured  along  the  right  abutment  at  El. 501*5  (in  feet)  near  the  toe  of 
the  dam  (Page  36). 

Observations  and  Dimension  of  the  Path  for  a  Second  Leak:  A  second 
leak  was  observed  between  10:00  a.m.  and  10:30  a.m.  at  El. 5195,  approxi¬ 
mately  1»6  meters  (150  ft)  above  the  first  leak  at  E1.50l»5.  By  10;l*0  a.m., 
an  eyewitness  observed  that  the  leak  had  developed  into  an  erosion  tunnel, 

about  1.8  meters  (6  feet)  in  diameter,  extending  into  the  embankment 
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(Page  36).  Unfortunately,  that  was  the  only  quantitative  description  of 
che  erosion  channel  provided  in  the  report.  Moreover,  available  photo¬ 
graphs  of  the  eroding  dam  for  different  times  do  not  provide  a  useful 
scale  that  could  be  used  to  quantify  the  dimensions  of  the  erosion 
channel. 

Progression  of  the  Erosion  Channel:  Seed  and  his  co-vorkers 
(U.S.  Department  of  the  Interior,  1976  Appendix  D)  described  the  pro¬ 
gression  of  flow  and  the  development  of  the  erosion  channel  through  the 
Teton  dam  embankment.  This  progression  is  illustrated  in  Figures 
FI  and  F2.  They  concluded  that  flow  first  occured  along  the  base  of  the 
dam,  resulting  in  the  first  leak  observed  at  El. 50U5  ft.  This  flow 
path  later  shortened  and  allowed  water  to  flow  out  of  the  embankment  at 
El. 5195,  resulting  in  a  second  leak.  This  second  leak  enlarged,  into 
the  1.8  meter  (6  ft)  in  diameter  erosion  channel  dencribed  above,  by 
10: 1*0  a.m.  on  June  5,  1976. 

It  may  be  noted  that  the  maximum  reservoir  level  on  the  eve  of  the 
failure  was  given  at  El.  $300  ft.  Thus  the  head  of  water  above  the  erosion 
channel  at  El. 5195  ft.  was  approximately  32  m  (105  ft). 

Calculations  to  Estimate  Erosion  Rates  from  Field  Observations: 
Computations  were  performed  to  obtain  estimates  of  the  erosion  rate  from 
the  information  given  above.  These  computations  used  the  two  procedures: 
l)  A  procedure  similar  to  the  one  outlined  in  the  text  of  this  report 
using  assumptions  for  the  weight  of  material  eroded,  estimates  of  the 
duration  of  erosion,  and  assumptions  for  the  erosional  surface  area;  2) 

A  procedure  outlined  by  Perry  (1982)  using  a  crack  recession  parameter, 

6,  defined  as  the  maximum  erosion  rate  tt P,v  divided  by  the  dry  density 
of  the  core  material.  This  recession  parameter,  6,  simply  describes 
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the  linear  enlargement  of  an  erosion  channel.  Details  for  using  these  two 
procedures  to  calculate  the  quantity  of  material  eroded  from  the  Teton 
Dam  are  outlined  below. 

Calculation  of  the  Erosion  Rate  as  Outlined  in  this  Report  (Method  A): 
This  procedure  requires  knowledge  of  the  volume  of  material  eroded  from 
the  core,  the  weight  of  this  volume,  the  time  for  this  erosion  to  occur, 
and  the  erosional  surface  area.  For  these  computations  it  was  assumed 
that  l)  the  length  of  the  erosion  channel  was  90  meters  measured  along 
the  base  of  the  core  as  shown  in  Fig.  F2d,  2)  the  shape  of  the  erosion 
channel  was  circular  with  a  diameter  of  1.8  meters,  and  3)  the  duration 
of  erosion  ranged  from  a  maximum  time  of  820  minutes  (9:00  p.m.  on  June 
1*  to  10:1*0  a.m.  on  June  5,  1976)  to  a  minimum  time  of  175  minutes  (7:1*5 
a.m.  to  10:1*0  a.m.  on  June  5,  1976). 

Volume  of  Material  Eroded  from  the  Core  by  10;  1*0  a.m. 

Vc  «  Ac  Lg  «  229  m* 

where:  Vc  ■  Volume  of  material  eroded  from  the 

core  at  10:1*0  a.m.; 

Ac  *  Cross-sectional  area  of  erosion  tunnel 
observed  at  10;1*0  a.m.  (equal  to  i*  (l.8)2 
or  2.51*  sq.  meters);  T* 

Lc  ■  Approximate  length  of  the  core  from  upstream 
to  downstream  at  El. 5195  (equal  to  90  meters). 

Weight  of  Core  Material  Lost  by  10: 1*0  a.m. 

Wc  ■  Pc  Vc  ■  362  Mg 

where:  Wc  ■  Total  weight  of  core  material  lost  by 

1C;!*0  a.m. ; 

p  *  Mean  dry  density  obtained  for  the  core 
material  (tone  l)  during  construction 
(equal  to  1.583  H|). 


.  „ j  -n  A^oded.  from  the  core  by 
V  =  Volume  of  material  eroaea  ir® 

C  i0:40  (equal  to  225  m  >  • 


fipt.-i mate  of  Rate  of 


uoight.  Erosion  Until  10:**0  a,B‘ 


wc  i  p  i  H?. 


0.1+3  <  P  <  1.98  Mg. 

min 


where : 


p  =  Rate  of  weight  erosion; 

vr  -  Total  weight  of  core  material 
C  lost  (equal  to  356  Mg); 

Tmax*  Elapsed  time  between  9:00  p.m.,  on 

June  4,  and  10:40  a.m.,  on  June  5,  197 

(equal  to  840  minutes); 

T  +  Elapsed  time  between  7:45 

T“in  a.m!,  June  5,  1976  (equal  to  175  minutes). 


Frrn*""  "»tp  at  10:U0  a,a- 


Pmin  5  fec  5 


> . 08  <  cc  <  0.39  g/cm2-  min 


where : 


t  =  Estimate  of  erosion  rate  at  10:40  a.m. 
(g/cm2-  min); 

Pmln.  Minimum  «Ummt.d  r.t.  of  nUM  «ro.ioo 

(0.43  Mg/min); 

P  ,  Maximum  estimated  rate  of  weight  erosion 
*“  (1.98  Mg/min); 

a  -  Erosional  surface  area  !°di«eter 

meters  through  the  core  and  with  a  diameter 

of  1.8  meters  at  10:40  a.m.  (equal  to  5-1 
•  10*  cm1). 
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It  is  interesting  to  note  that  this  calculated  range  of  erosion 
rate  values  (0.08  to  0,39  g/cm2-  min)  is  nearly  identical  to  maximum 
erosion  rate  values  obtained  from  Laboratory  Triaxial  Erosion  Tests  per¬ 
formed  on  Teton  Dam  specimens  (iteterial  E).  These  measured  laboratory 
erosion  rates  ranged  between  0.13  and  0.35  g/cm2-  min. 

Calculation  of  Erosion  Rate  as  Outlined  by  Perry  (1982)  (Method  B): 
The  maximum  erosion  rate  values  obtained  from  Laboratory  Triaxial  Erosion 
Tests  presented  in  this  report  may  be  used  to  compute  the  recession  of 
the  wall  of  a  crack  due  to  erosion  using  the  following  expression 
(Perry  1982) : 


6  =  —p —  ft /day) 

cm/min 


(FI) 


where: 


6  *  recession  of  wall  of  crack  due  to 
erosion,  ft/day; 


^max  =  maximum  erosion  rate,  g/cm2-  min;  and 
p  *  dry  density  of  soil,  g/cm3. 


This  relationship  between  the  crack  recession  rate,  6,  and  the  maximum 
erosion  rate,  cm«y.  is  plotted  jot  Figure  F3  for  various  values  of  dry 
density.  The  recession  rate  calculated  using  Equation  FI  may  represent 
an  upper  bound  for  erosion  for  a  dam  having  a  downstream  filter  because 
it  does  not  take  into  account  accretion  of  eroded  material  at  the  inter¬ 
face  between  the  cracked  core  and  the  downstream  filter.  On  the  other 
hand,  this  recession  rate  may  represent  a  lower  bound  for  erosion  of  a 
crack  that  is  enlarging  due  to  both  wall  erosion  and  the  collapse  of 
the  embankment  into  the  erosion  channel.  Thus,  unless  corrections  are 
applied  for  filters  and  the  rate  of  embankment  collapse,  Equation  FI 
will  accurately  model  stable  crack  growth  due  to  wall  erosion  in  an  em¬ 
bankment  unprotected  by  a  downstream  filter. 
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Reviewing  the  design  drawings  and  the  field  performance  during 
the  failure,  it  is  clear  that  the  Teton  Dam  embankment  was  not  adequately 
protected  by  a  downstream  filter.  Further,  from  laboratory  tests  and 
from  observations  made  during  the  failure  of  the  dam,  it  is  clear  that 
the  core  material  collapsed  into  the  erosion  channel.  Thus,  Equation  FI 
will  appear  to  give  a  lower  bound  for  values  of  the  recession  of  the  wall 
of  a  crack  in  the  Teton  Dam. 

Example  Calculation  of  Maximum  Erosion  Rate 

The  two  methods  presented  above  may  be  used  to  back  calculate  the 
erosion  rate  that  took  place  during  the  Teton  Dam  failure. 

First,  the  recession  rate,  6,  may  be  computed  for  the  Teton  Dam 
assuming  that  the  crack  enlarged  from  a  negligible  thickness  and  only 
in  the  vertical  direction.  Further,  it  is  reasonable  to  assume  based 
on  failure  observations  that  the  crack  had  a  constant  width  of  1.8  meters 
(6  ft)  and  that  it  increased. only  in  thickness  from  0  meters  to  1.8 
meters  (0  ft.  to  6  ft.)  during  failure.  This  implies  that  the  crack 
eroded  as  a  constant  width  rectangle.  This  is  a  logical  assumption 
if  the  crack  was  very  wide  compared  to  it's  thickness.  The  time  for 
recession  of  the  crack  ranged  from  a  maximum  of  0.57  days  (820  min)  to 
a  minimum  of  0.12  days  (175  min).  Thus  the  recession  rate,  6,  ranged 
from  3.2  m/day  (10.35  ft/day)  to  15  m/day  (50  ft/day). 

Second,  it  is  necessary  to  know  the  average  dry  density,  p,  of 
the  eroding  dam  material.  From  field  tests  it  is  known  that  the  average 
dry  density  of  the  core  material  (zone  l)  was  1.583  Mg/m*. 
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Range;  of  Estimated  Maximum  Erosion  Rates:  Using  these  values  of 
the  recession  rate  and  the  core  material  dry  density,  a  range  of  values 
of  the  maximum  erosion  rate,  em av ,  can  be  read  off  from  Figure  F3  (Method 
B).  The  resulting  values  for  the  maximum  erosion  rate  is  from  0.36  to 
1.7  g/cm2-  min  which  is  about  5  times  greater  than  the  range  of  values 
computed  using  the  data  and  procedures  developed  in  this  report  (Method 
A).  This  wide  spread  between  calculated  values  using  Method  B  and  the 
values  computed  using  Method  A  may  be  due  to  the  effect  of  the  implicit 
assumption  made  for  the  values  obtained  using  Method  B  that  the  erosion 
channel  enlarges  only  in  one  direction.  Thus,  Method  B  should  be  applied 
to  estimate  the  erosion  rates  in  dams  where  the  crack  enlarges  such  that 
the  height  and  width  of  the  crack  remain  approximately  the  same  as  in 
the  case  of  circular  erosion  channel. 

Design  Use 

These  computations  indicate  that  values  for  the  maximum  erosion 
rate  as  measured  in  the  Triaxial  Erosion  Test  are  reasonable  and  may 
be  used  in  design.  However,  the  procedure  for  estimating  the  erosion 
rate  of  a  eroding  dam  outlined  in  this  Appendix  requires  that  a  variation 
of  crack  shape  with  time  must  be  assumed.  Method  A,  outlined  in  this 
Appendix,  will  be  most  useful  for  a  wide  variety  of  crack  shapes.  Method 
B,  however,  will  be  most  useful  if  a  crack  is  assumed  (or  measured)  to 

i 

be  enlarging  faster  in  one  direction  relative  to  it's  other  direction. 

An  example  of  such  a  case  is  a  long  crack  at  the  interface  between  the 

, 

core  and  bedrock,  or  a  long  crack  between  compaction  lifts.  For  such 
crack  geometries,  the  recession  rate  analysis  (Method  B)  will  simplify 
computations  and  should  result  in  reasonable  values  for  estimated  crack 
enlargement  with  time. 
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Recommendations 


The  autho  o  believe  that  analyses  of  the  rate  of  dam  core  erosion 
could  become  mor  credible  if  better  eyewitness  reports  of  embankment 
dam  core  erosion  are  provided  in  the  literature.  As  a  minimum  these 
reports  should  in,  lude  information  on:  l)  Fluid  flow  rate  values;  2) 

An  estimate  of  the  flow  channel  shape  or  span  length;  3)  The  weight  of 
solids  being  eroded  as  a  function  of  time;  and  k)  Estimates  of  the  en¬ 
largement  of  the  erosion  channel  with  time. 
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APPENDIX  G:  NOTATION 


SYMBOL 
As(  t) 

As(0) 

Ax(t) 

Ax(0) 

d 

h 

hb 

i 

L 

1e 

P 

Q 


description 

Eroded  surface  area 
as  a  function  of  time 

Eroded  surface  area 
at  time  equal  zero 

Cross  sectional  area 
of  crack  in  specimen 
as  a  function  of  time 

Cross  sectional  area 
of  crack  in  specimen 
at  time  equal  zero 

Hydraulic  diameter 

Thickness  of  crack¬ 
forming  blade  (crack 
thickness) 

Head  of  water  applied  to 
the  top  of  the  sample 

Head  of  water  applied  at 
the  end  of  outflow  line 
(equal  to  zero  in  these 
tests) 

Hydraulic  gradient 

Initial  length  of  the 
specimen 

Eroded  (crack)  length 
Rate  of  weight  erosion 
Plow  rate 


UNITS 

2 

cm 

cm2 

cm2 


cm 

cm 

m 

m 

m/m 

cm 

cm 

g/min 
cm* /min 


R 


e 


S 


Reynold’s  number 

Ratio  of  solids  eroded 
to  volume  of  flow  collected 
at  any  time  t 
(Parts  per  million) 


mg/liter 

(ppm) 


G1 


(*>  • 


SYMBOL 

T 

t 

Wi 


w 


OC 


V 

p 

a3c 

T 

* 


DESCRIPTION 

UNITS 

Elapsed  time 

minutes 

time 

minutes 

Dry  weight  of  specimen 
at  start  of  test 

grams 

Dry  weight  of  specimen 
after  erosion  test 

grams 

Width  of  slot-forming 
blade  (crack) 

cm 

Percent  erosion  per  unit 
time 

J/min 

Erosion  rate 

g/(cm2-min) 

Viscosity  of  Water 

N-Sec 

m2 

Dry  density  of  specimen 

g/cm* 

Triaxial  Cell  Pressure 

KN/m* 

Fluid  Shear  Stress 

N/m2 

Flow  parameter 

— 
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